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ABSTRACT 

We present a new library of fully-radiative shock models calculated with the MAPPINGS in shock 
and photoionization code. The library consists of grids of models with shock velocities in the range 
v s = 100-1000 km s' 1 and magnetic parameters B/^/n of 10~ 4 -10 ^G cm 3 / 2 for five different atomic 
abundance sets, and for a pre-shock density of 1.0 cm' 3 . Additionally, Solar abundance model grids 
have been calculated for densities of 0.01, 0.1, 10, 100, and 1000 cm' 3 with the same range in v s and 
B/i/n. Each model includes components of both the radiative shock and its photoionized precursor, 
ionized by the EUV and soft X-ray radiation generated in the radiative gas. We present the details 
of the ionization structure, the column densities, and the luminosities of the shock and its precursor. 
Emission line ratio predictions are separately given for the shock and its precursor as well as for the 
composite shock+precursor structure to facilitate comparison with observations in cases where the 
shock and its precursor are not resolved. Emission line ratio grids for shock and shock+precursor are 
presented on standard line ratio diagnostic diagrams, and we compare these grids to observations of 
radio galaxies and a sample of AGN and star forming galaxies from the Sloan Digital Sky Survey. 
This library is available online, along with a suite of tools to enable the analysis of the shocks and the 
easy creation of emission line ratio diagnostic diagrams. These models represent a significant increase 
in parameter space coverage over previously available models, and therefore provide a unique tool in 
the diagnosis of emission by shocks. 

Subject headings: hydrodynamics - shock waves - ISM: abundances,- Galaxies: Nuclei, Galaxies: 
Seyfert - infrared: ISM, Ultraviolet: ISM, X-rays: ISM 



1. INTRODUCTION 

Supersonic motions are a common phenomenon in our 
complex and fascinating Universe. The kinetic energy of 
such motions will almost inevitably be eventually dissi- 
pated through radiative shocks. Cloud-cloud collisions, 
the expansion of H II regions into the surrounding inter- 
stellar medium, outflows from young stellar objects, su- 
pernova blast waves, outflows from active and starburst 
galaxies, and collisions between galaxies are all examples 
of astrophysical situations in which radiative shock waves 
provide an important component of the total energy bud- 
get and may determine the line emission spectrum. 

In this paper we do not consider very slow (molecu- 
lar) shocks or indeed, the faster atomic shocks for which 
the theory ha s b een well-d e velop ed by such authors as 
iDopital (|1977t ) or iRavmoridl (|1979j ). Rather, we consider 
only fast shocks, where the ionizing radiation generated 
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by the cooling of hot gas behind the shock front gener- 
ates a strong radiation field of extreme ultraviolet and 
soft X-ray photons, which leads to significant photoion- 
izing effects. The detailed theory of steady-flow pho- 
toionizing (or au to-ionizing ) shocks was developed by 
Sutherland et alJ (|1993f ) and lDopita fc Sutherland i|1995l . 
19961 ) (hereafter DS95 and DS96). A detailed text- 
book development of the t heory of shocks is given in 
iDopita fc Sutherland! (|2003D . 

In photoionizing shocks, the flux of the ionizing ra- 
diation emitted by the shock increases in approximate 
proportion to the energy flux through the shock (oc v 3 ). 
The ratio of this flux in advance of the shock to the 
pre-shock density, classified as the ionization parameter, 
determines the velocity of the photoionization front that 
is driven into the pre-shock gas. At low values of the 
ionization parameter, the velocity of the photoionization 
front is lower than the velocity of the shock and the ion- 
izing photons are absorbed in the immediate vicinity of 
the shock front. The effect of this is to change the ion- 
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Abundance Sets (by number wrt Hydrogen) 
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Fig. 1. — Ionizing spectra generated by v s =250, 500, 750 and 
1000 km s" 1 shock models with n=1.0 cm" 3 and solar abundance. 



ization state of the gas feeding across the shock front. 
As the velocity of the shock increases, the emitted ioniz- 
ing flux, and therefore the velocity of the photoionization 
front, increases rapidly. At shock velocities above a cer- 
tain limit (v a « 170 km s _1 ), the ionization front veloc- 
ity exceeds that of the shock (and is supersonic with re- 
spect to the pre-shock gas) and the photoionization front 
detaches from the shock front as an R-Type ionization 
front 1 . This front expands to form a precursor H II region 
ahead of the shock. At the highest shock velocities, the 
photoionized precursor emission may come to dominate 
the optical emission of the shock, and the global radia- 
tive shock spectrum provides a rich mixture of emission 

R-Type ionization fronts are defined as being supersonic with 
respect to the gas ahead of the front, (subsonic ionizati on fronts 
are classified as D-Type) sec McKec & Hollcnbach (1980) 
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Fig. 2. — Ionizing spectra of d s =400 and 1000 km s -1 shock 
models shown on vF v scale. 

lines from both high- and low-ionization species. 

Appreciable photoionization effects are also produced 
in the gas behind the shock front, near to the recombi- 
nation region of the shock. However, in this shock region 
the velocity of the ionization front is much slower due 
to the gas compression through the shock front and the 
ionization front stays trapped in the recombination re- 
gion to much higher shock velocities. When magnetic 
effects are negligible, the compression of the gas in the 
shocked region is proportional to the square of the Mach 
number in the pre-shock gas, M = v s /cn, where c\\ is 
the sound speed in the pre-shock gas. When the mag- 
netic field pressure in the post-shock gas dominates over 
the gas pressure, the compression factor is determined by 
the Alfven Mach Number, Ma, which is the ratio of the 
shock velocity to the Alfven velocity; v\ = (B 2 /47rp) 1 / 2 , 
where B is the transverse component of the pre-shock 
magnetic field, and p is the pre-shock mass density. Thus 
the presence of magnetic fields act to limit the compres- 
sion through the shock. In order to account for this effect 
in computations of the shock structure, DS96 developed 
the concept of a magnetic parameter; T^/^fn., where n is 
the pre-shock particle number density. 

The grid of low density photoionizing shock models 
described in DS96 has proved to be a valuable resource 
to the astronomical community for assessing the role 
of photoionizing shocks in a range of astrophysical ob- 
jects. This utility is increased by their availability in 
elec tronic form from the AAS CD-ROM Series, Vol 7. 
(see iLeitherer et ailll996f) . These models have most of- 
ten been applied in studies of the narrow line regions 
(NLR) of active galaxies, developing the debate on the 
relative roles of shock and/or photoionization excitation 
of the NLR, and the question of the radiative versus me- 
chanical energy output of active galactic nuclei (AGN). 
Other uses of the models have includ ed studies of ion - 
ized gas around high velocity clouds (|Fox et alj [2004) . 
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Fig. 3. — The hydrogen ionization structure, temperature profile and density profile of a shock and precursor. The vertical axis represents 
i) the ionization fraction X; on (values < 1) for the ionization structure curves, ii) the electron temperature profile (T e ) in degrees Kelvin, 
shown as the upper line plotted in the figure, and iii) the hydrogen density (nu) in cm -3 , shown by the dotted line in the shock region 
(right) panel of the figure. As npj is constant in the precursor (left) it is not plotted. The horizontal axis represents the time since passage 
of the shock front, with positive values for the shock structure shown in the right of the figure, and negative values for the precursor region 
where the shock is yet to arrive. The numerical labels indicated on the temperature profile of the shock correspond to different regions of 
the shock structure as described in the text. The model plotted here is our fiducial model with solar abundances, equipartition magnetic 
field, a precursor density of njj = 1 cm -3 , and shock velocity of i) a =500 km s -1 , as labelled in the upper right corner of the figure. See the 
electronic edition of the Journal for a color version of this figure. 

and radiation fields associated with shocks is necessary 
in order to be able to more closely model these com- 
plex physical systems, and to be able to draw deeper 
conclusions. Indeed the 1-D, steady-flow nature of these 
models is acknowledged as their greatest limitation, since 
all fast shocks are subject to thermal instability, and 
the generation of secondary shocks within the cooling 
zone. To address these issues, much effort is being put 
into higher dimensional numerical simulations of radia- 
tive shocks, in particular in the area of understanding 
the role of local thermal and dynamical instabilities on 
both the shock structur e and on the emergent spectrum 
([Sutherland et al.ll2003f ). However, given the complexi- 
ties inherent in their computation, and the large amount 
of supercomputer time required, complete grids of model 
predictions from full 3-D models are some time off, and 
until then, modeling will most likely favour simulations 
of individual and specific cases. 

Thus, the simpler 1-D models remain a useful tool 
for comparing with observations of both individual ob- 
jects and in the investigation of the mean parameters 
of groups of objects observed in surveys. Furthermore, 
these models serve as a stepping stone towards higher 
dimensional models. For example, under some simplify- 
ing assumptions the radiative properties of 1-D models 
may be mapped onto those of 3-D hydrodynamical mod- 
els. In particular, thermally unstable shocks of a given 
velocity, v s , behave si milarly to steady-flow sh ocks with 
shock velocities ~ |u s ([Sutherland et al.ll2003l ). Bearing 
in mind these effects thermal instability, we advise that 
care should be taken in the use of 1-D models to derive 
physical parameters in individual cases, or at least allow 
for this factor of 2/3. 
In this paper we present a new library of 1-D steady 



and the intergalactic medium (jShull et al. 2003). 

In DS95 these models were applied to optical line ratios 
observed in NLRs. They found that Seyferts have values 
close to the shock + precursor predictions and that Low 
Ionization Emission Line Region (LINER) galaxies fall 
within the range of the shock-only models. In I Allen et al.l 
( 1998) we highlighted the importance of UV line ratios 
which help separate shock and photoionization models, 
and we defined a set of UV line ratio diagnostic diagrams. 
These UV diagnostics wer e applied to individual objects 
M 87 (iDopita et all H997TI . NGC 5 728 and NGC 564 3 
(|Evans et al.lll999ft and NGC 1068 (|Groves et al.ll200l . 
NGC 2992 was also investigated in detail with optical 
diagn ostics and other multi- wavelength data (|Allen et al.l 
[19991) . 

The DS95 model grid has also been extensively 
used by ot her authors i n stud ies of AGN. iBest et all 
(|2000t ) and llnskip et all (|2002l ) used optical and UV 
line ratios to show that shocks associated with ra- 
dio sources in z ~ 1 3CR radio galaxies can dom- 
inate the kinematics and ionization as the radio 
source expands through the interstellar medium of 
the host galaxy. These models have also application 
to the analys is of emission lines in h i gh redshift ra- 
dio galaxies dVillar-Martin et alj Il997t iMaxfield et all 
2002t iDe Breuck et all 120001: IVillar-Martm et al.l l20Qg 



Iwamuro et al.l 120031 : iReuland et al.ll2007l ). and in stud- 



ies of indivi dual low redshift a ctive gala xies including 



NGC 4151 (iNelson et al I [200(3) . M 51 ([Bradley et al 
2001 . M 87 (ISabra et alJ 120031) . Mkn 78 dwhittle et a" 



200il . 3C 299 (iFeinstein et aj.lH^9l ). NGC 2110 and 



NGC 5929 (iFerruit et al.lll999D . 

Studies of individual objects make it obvious that 
higher dimensional modeling of the physical structures 
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flow photoionizing shock models, calculated with an 
updated version of the shock modeling code, MAP- 
PINGS in. The library is an improvement upon and ex- 
tension of the DS95 & DS96 models, and includes a range 
of chemical abundance sets, pre-shock densities from 0.01 
to 1000 cm -3 , velocities up to 1000 km s _1 , magnetic 
fields from 10~ 10 to 10 -4 G, and magnetic parameters 
(B/Vn) from 10~ 4 to 100 ^G cm 3 / 2 . These new model 
grids therefore supersede the previous models, and are 
designed to be of maximum utility to observers for com- 
paring observations with photoionizing shock models. 

As well as reiterating the technique used in MAP- 
PINGS to create the shock models, we discuss in detail 
the radiation fields generated by shocks, and the resulting 
ionization, density and temperature structure of shocks 
for the range of parameters considered. We also present 
the resulting model grids on a range of ultraviolet, op- 
tical and infrared emission-line diagnostic diagrams, dis- 
cussing the effects of the various parameters upon the 
grids, and demonstrate their usefulness and applicabil- 
ity by comparing these grids with emission-line galaxy 
sample from the Sloan Digital Sky Survey (SDSS). 

We introduce the MAPPINGS in online shock library, 
the main access point for the results of the shock models, 
which includes both the resulting emission lines and col- 
umn depths from the models, as well as the full results 
of the MAPPINGS in models for detailed analysis. Also 
included as part of the library are the analysis tools used 
within this work, providing a simple and easy mechanism 
to determine the shock diagnostic power of any emission- 
line ratios. 

2. MODELING TECHNIQUE 

The models presented here have been calculated with 
the MAPPINGS in shock and photoionization mod- 
eling code, version Illq. This is an updated ver- 
sion of the MAPPINGS II c ode which was described in 
ISutherland fe Dopital (fl993). The main enhancements of 
version III include the use of a higher resolution radia- 
tion vector, and a more stable scheme for choosing time- 
steps which allows computation of higher velocity shocks. 
Other improvements include the explicit inclusion of all 
the He II emission lines, and a corrected computation 
of the neutral oxygen emission. Further details on the 
use of and the chronology of the development of MAP- 
PINGS models may be found at the MAPPINGS Online 
web pages 2 . Various improvements to the calculation 
of photoion ization processes have be en implemented in 
version III (|Groves et al.ll2004al l2007h although the dust 
effects discussed there have not been considered in these 
shock models, which are fundamentally dust-frce. This 
approximation is probably a physical one, since in the 
faster shocks dust will be effectively destroyed by grain- 
grain collisions, through both shattering and spallation, 
and by thermal sputtering. The ph ysics of these pro- 
cesses h ave been discussed in det ai l iDraine fc Salpeterl 
1979h : iDwek Foster fc Vancural (119961): Pones et all 



1994 I1996D; iPineau de Forets fc Flower! (|1997f) and 
iPerna fc LazzatiT i 20021 ). We have in any case run models 



with depleted abundances to facilitate comparison with 
the earlier DS96 high- velocity shock modeling, so our dif- 
ferent model grids can be used to provide some estimate 

2 http: //www. if a.hawaii . edu/~ {}kewley/Mappings 



of the effect of dust depletion on the output spectra. 

MAPPINGS in uses exactly the same computational 
recipe for calculating the flow solution of these one- 
dimensional radiative shock models as described in detail 
in DS96 for MAPPINGS II. Briefly, the hydrodynamics 
of the flow are derived from the Rankine-Hugoniot jump 
conditions including the magnetic field terms. The four 
Rankine-Hugoniot conditions, equations [T] - [4j represent 
the conservation of mass, conservation of momentum, the 
condition that the magnetic field is locked i nto the ion- 
ized p lasma, and conservation of energy (see iHeng et all 
(2007) for a generalization of these conditions). These 
equations connect two points in the flow at times t and 
to with density, velocity, pressure, internal energy, and 
transverse magnetic field p, V,P, [/, B, and po, Vq,Pq, 
Uo, Bo, respectively. A is the mean cooling rate of the 
plasma over the time step. As set out in DS96, by choos- 
ing the subscripted flow variables as the point immedi- 
ately in front of the shock, the equations can be reduced 
to a quadratic expression for the flow velocity. By choos- 
ing appropriate fractions of the geometric mean of the 
plasma timescales, the physical conditions in the plasma 
can be followed smoothly. 
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The main calculation to be performed at each time 
step in these shock models is that of the ionization bal- 
ance and the corresponding mean cooling rate, A. The 
cooling and radiative emission is calculated using a very 
large atomic data set which allows treatment of all ion- 
ization stages of cosmically abundant elements up to fully 
ionized nickel. Using this atomic data set, the rate equa- 
tions for non-equilibrium ionization, recombination, ex- 
citation, and radiative transfer and cooling are solved at 
each time step of the flow. Details of the various colli- 
sional and radi ative processes included in MAPPINGS 
are described in ISutherland et all (|1993l ). 

Modeling of all fast radiative shocks is necessarily an 
iterative process because the detailed structure of the 
shock depends on the ionization state of the precursor 
gas entering the shock front. To calculate a fully self- 
consistent model of a radiative shock a number of it- 
erations are required. Firstly an initial shock model is 
calculated using an estimate of the ionization state of 
the precursor. Then the photoionized precursor is cal- 
culated using the ionizing radiation field generated by 
the shock. This process is then repeated, updating the 
ionization state of the precursor gas at each iteration. 
The models presented here employ four such iterations, 
which was found to be sufficient to allow the temperature 
and ionization state of the precursor gas to stabilize at a 
constant value. 

In the final iteration of each model, the computation 
of the shock structure is allowed to proceed until the gas 
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Fig. 4. — The hydrogen ionization structure, temperature profile and density profile of a shock and precursor for models with shock 
vcocities of v B = 200 and 1000 km s — 1 and equipartition magnetic field. The axes are as described in Figure [3] See the electronic edition 
of the Journal for a color version of this figure. 

has cooled to 1000K, below which no further significant 
emission in the considered species is produced. The com- 
putation of the precursor structure is terminated when 
the ionized fraction of hydrogen falls to below 1 percent. 

Precursor components were computed for each individ- 
ual model, in contrast to DS96 who computed precursors 
only for each value of the shock velocity. The precursor 
ionization does not generally depend on the magnetic 
parameter but here we chose to compute each individ- 
ual precursor to ensure full self consistency between the 
shock and its precursor. 

2.1. Model Grid Input Parameters 

Each individual model in our shock library is defined 
by five physical parameters; the pre-shock density, n, 
the shock velocity, v s , the pre-shock transverse magnetic 
field, B and the set of atomic abundances. The library 
itself is organized into two main groups of model grids: 

• First, complete grids of models calculated for five 
different atomic abundance sets; (depleted) Solar 
and 2 x Solar as used by DS96 and w h ich ar e based 
upon the older lAnders fc Grevessel (|1989f ) abun- 
dance set; a Solar abundance set based upon the 



lAsplund. Grevesse fc Sauval 1 2005f ) abundances as 
listed in iDopita et al! I (|2005lT which is referred 
to here as the 'dopita2005' abundances; and an 
LMC and an SMC ab undance set as given by 
iRussell & Dopital (|1992f L The abundances of the 
individual elements (by number with respect to hy- 
drogen) are listed for each of these abundance sets 
in Table [T] Each of these model grids assume a 
fixed pre-shock density of n=l cm -3 , and consists 
of a set of models with shock velocities covering 
the range 100 up to 1000 km s _1 in steps of 25 km 
s^ 1 , and magnetic fields of 10~ 4 , 0.5, 1.0, 2.0, 3.23, 
4.0, 5.0 and 10.0 /iG. As the pre-shock density is 
unity, the corresponding magnetic parameters are 
B/Vn= 10~ 4 , 0.5, 1.0, 2.0, 3.23, 4.0, 5.0 and 10.0 
/xG cm 3 / 2 . 

• Second, grids of models with solar abundance (as 
used by DS96), calculated for densities of 0.01, 0.1, 
1.0, 10, 100 and 1000 cm" 3 . As with the first set, 
each of these model grids covers shock velocities of 
100 up to 1000 km s _1 in steps of 25 km s" 1 , and 
have the same magnetic parameters of 10~ 4 , 0.5, 
1.0, 2.0, 3.23, 4.0, 5.0 and 10.0 fiG cm 3 / 2 . How- 
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Fig. 5. — The hydrogen ionization structure, temperature profile and density profile of models with different magnetic parameter. Models 
of v s = 500 km s _1 are shown for the magnetic parameters of 0.0001 and 10 fiG cm 3 / 2 . The axes are the same as for Figure [3] See the 
electronic edition of the Journal for a color version of this figure. 



ever, as the magnetic field values required to obtain 
these magnetic parameters are of course different 
for each density, additional models were computed 
in order to be able to also compare models of dif- 
ferent densities with the same transverse magnetic 
field. These are B~ 1CT 3 , ~ 1CT 2 , ~ 10 -1 , 1.0, 10 
and 100 /iG, calculated for each density. 

The magnetic field, B, and magnetic parameter, H/y/si, 
values are chosen so as to cover the extremes expected 
in the ISM, while also sampling more finely the mag- 
netic field strengths which are near equipartition. Under 
equipartition conditions the magnetic pressure is equal 
to the thermal pressure, and the Alfven speed is approx- 
imately equal to the gas sound speed. Pressure equiparti- 
tion occurs for Bq/47t ~ n$kT$ where Bo is the transverse 
magnetic field, and no, Tp are the pre-shock densities and 
temperatures (jDopita fc Sutherland! f2Q03h . This condi- 
tion is satisfied for magnetic parameters B/y 7 ?! ~ 3 — 5. 
The value B/y/n ~ 3.23 was chosen as the nominal 
equipartition value. 

The input parameters for the all of the models are pro- 
vided in Table [3] Each row of the table represents a ve- 
locity sequence of models at a given abundance, density 



and magnetic field. The table is organized into ten sets 
of models. The first five sets of models represent the 
model grids for the five different abundances, and the 
following 5 sets of models are the solar abundance mod- 
els for different input densities and magnetic fields. The 
number and range of model input parameters allows the 
construction of various 2-D, or higher dimensional, grids 
of models. 

3. IONIZING RADIATION GENERATED BY THE SHOCK 

The ionizing radiation produced in the cooling zone 
behind the shock front shock is mostly composed of 
thermal bremsstrahlung (free-free) continuum and res- 
onance lines arising from many different elements and 
ionic stages. The underlying exponential shape of the 
continuum is emphasized in Figure [T] where we show 
the ionizing spectra generated in the n=1.0 cm -3 solar 
abundance shock models. As can be seen, higher ve- 
locity shocks result in harder and more luminous ioniz- 
ing spectra, with the spectral slope in the log-linear plot 
clearly flattening with increasing velocity. Note that the 
histogram nature of the figure also reveals the high en- 
ergy (frequency) resolution of the MAPPINGS spectral 
vector. Figure [2] shows the ionizing spectra of 400 and 
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Fig. 6. — The hydrogen ionization structure, temperature profile and density profile of models with different pre-shock density. Models 
of v a = 500 km s — 1 and equipartition magnetic field are shown for pre-shock densities of 0.1, 10 and 1000 cm -3 . The axes are the same as 
for Figure [3] See the electronic edition of the Journal for a color version of this figure. 

1000 km s -1 shocks on a vF v scale, illustrating both 
the stronger UV fluxes and harder X-ray spectra gen- 
erated in high velocity shocks. This figure also shows 
a prominent low-temperature bound-free continuum of 
hydrogen, produced in the cool, partially-ionized zone 
of the recombination region of the shock, and the strong 
hydrogen two-photon continuum produced mostly by the 
down-conversion of Lya photons trapped in this same re- 
gion of the shock structure. Also present, though to a 
much weaker scale, is the bound-free continuum arising 
from the heavier elements, with the helium continuum 
the most obvious. 
The strength of the ionizing field is a strong function 



of the shock velocity, but does not significantly depend 
on either the atomic abundance or the magnetic field. 
This is because the ionizing field is dominated by the 
bremsstrahlung radiation, whose strength is controlled 
by the temperature and density of the radiative zone, 
which is determined by the shock velocity and pre-shock 
density. As magnetic field support in this cooling zone is 
negligible, it has little affect on the emission, and, as hy- 
drogen dominates both electron and ion numbers for both 
the bremsstrahlung and free-bound emission, changes in 
metallicity have only a small impact. 

The hydrogen ionizing radiation flux, -LtjVi integrated 
for all energies hv > 13.6 eV, and over 2ir sr is listed 



Allen et al. 



§ 10 2 



Std2.abn V=0500 n=1.0 B=3.22640 



-60 -40 
TIME (10" s) 



2 3 
TIME (10" s) 



I- 10* - 



§ 10 2 



W 10 L 

V 




doplto2005.obn V=0500 n=l.0 B = 3.2264 



-140 -120 -100 - 80 - 60 -40 - 20 
TIME (10" s) 



10 20 

TIME (10" s) 



§ 10 2 



X 10-* - - 



Te„ 


smc.abn V=0500 n=1.0 B = 3.22640 




Te 


H II 




n 


H 1 


H II 










.... i .... i ... . 





-100 

TIME (10" s) 



10 15 
TIME (10" s) 



Fig. 7. — The hydrogen ionization structure, temperature profile and density profile of models with different atomic abundances. Models 
of v s = 500 km s~ 1 and equipartition magnetic field are shown for the 2 X solar, dopita2005 and SMC abundance sets. The axes are the 
same as for Figure [3] See the electronic edition of the Journal for a color version of this figure. 



in Table O It is found to scale almost exactly as the 
available enthalpy as: 



L uv = 2.44x10" 



vs 



100 km s' 



3.02 



ax " > 



24< hv < 54 eV band, Q(Hel), and the He II-ionizing 
band hv > 54 eV, Q(HeII) as a function of shock ve- 
locity. We also provide the corresponding R— Type ion- 
ization front velocity Vi on and the equilibrium electron 
ergs cm s temperatures, T e just ahead of the shock front. 



( 5 ) 

The ionizing fields may also be characterized in terms 
of their ionization parameter in the pre-shock gas. As 
in DS96 we give ionization parameters, Q, defined as 
the mean number of photons passing through unit area 
divided by the total pre-shock particle density; Q — 
Nphotons/nr (cm s _1 ). Table [2] lists Q in the H I- 
ionizing 13. 6< hv < 24 eV band, Q(H), the He I-ionizing 



4. SHOCK AND PRECURSOR STRUCTURES 

We now consider the ionization structures and the 
physical scales of both the shock and precursor compo- 
nents. Figure [3] shows the ionized and neutral hydro- 
gen structures (as labelled) of the shock and precursor 
components of the 500 km s _1 solar abundance model, 
along with the hydrogen density and electron tempera- 
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Fig. 8. — The helium ionization structure, temperature profile and density profile of models with solar abundance, equipartition magnetic 
field for shock velocities of v s = 200, 500, and 1000 km s . The axes are the same as for Figure[3] See the electronic edition of the Journal 
for a color version of this figure. 



ture profiles. Note that the hydrogen density is not plot- 
ted within the precursor region as it is set at the constant 
value of nn = 1.0 cm~ 3 . The profiles are plotted in the 
frame of the shock front, considered with respect to the 
time axis. This axis represents time since (or prior to) 
the passage of the shock front, and can be converted into 
a physical distance from the shock front via the shock 
velocity. 

In Figure [3] (and subsequent ionization structure Fig- 
ures m-[12]) , the structure of the precursor is shown in the 
left panel and the shock structure is shown in the right 
panel. The shock front is located at time t = 0. The 
time taken to reach equipartition in the electrons, ions, 



and un-ionized gas is assumed to be negligible, so that 
the shock is therefore unresolved and the density and 
electron temperature jump discontinuously at the shock 
front. Note also that the time axes are shown in units of 
10 12 s, but the scales and ranges are necessarily different 
for the shock and precursor panels of these diagrams, as 
well as between the diagrams with different parameters. 

As described in DS96 the a number of zones can be 
identified in the shock structures. These are indicated in 
Figure [3] and are, respectively: 

1. The ionization region. This region, just after the 
shock front, is in an ionization state below that 
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Fig. 9. — The carbon ionization structure, temperature profile and density profile of models with solar abundance, equipartition magnetic 
field for shock velocities of v s = 200, 500, and 1000 km s . The axes are the same as for Figure[3] See the electronic edition of the Journal 
for a color version of this figure. 



appropriate for the electron temperature assuming 
collisional ionization equilibrium (CIE), due to dis- 
continuous temperature jump at the front. After 
the passing of the front, the gas rapidly (barely 
resolved in the figures) adjusts from the pre-shock 
ionization state to the appropriate CIE state for the 
post-shock temperature. The rate at which it ad- 
justs depends upon the collisional ionization rates 
of these species. As discussed in DS96, this region 
has strong line emission due to this state of high 
temperature and relatively low-ionization. 

2. The high temperature radiative zone. This is the 



zone in which most of the EUV and soft X-ray flux 
is emitted. The ionization state is approximately 
in coronal equilibrium for the temperature of the 
region. Even though the cooling rate of this region 
is low due to the high temperature, it is the domi- 
nant contributor to the radiation field of the shock, 
assisted by the high-ionization, optically-thin state 
of the gas. 

3. The non- equilibrium cooling zone. Once the tem- 
perature of the post shock gas drops below ~ 10 6 
K, the cooling rate becomes very high and the re- 
combination timescales of a large number of ions 
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Fig. 10. — The nitrogen ionization structure, temperature profile and density profile of models with solar abundance, equipartition 
magnetic field for shock velocities of v s = 200, 500, and 1000 km s _1 . The axes are the same as for Figure|3] See the electronic edition of 
the Journal for a color version of this figure. 



becomes longer than the local cooling timescale. As 
a result, the ionization state lags and the plasma is 
in a higher degree of ionization than collisional ion- 
ization equilibrium would suggest. The relatively 
high ionization state for the temperature means 
that the collisional line emission is weak for these 
intermediate ionization species, and that the gas 
remains optically thin to the diffuse radiation of 
the previous region. 

4. The super-cooled zone is the region where photoion- 
ization starts to become important in determining 
the ionization balance. Initially, however, the ion- 



ization state is still too high to efficiently absorb 
the ionizing photons. This leads to an over-cooling, 
and the temperature falls below the value given by 
photoionization equilibrium. The width of this re- 
gion is dependent upon both the recombination and 
photoionization timescales of the ions. 

The photoab sorption and recombination zone. This 
is essentially an equilibrium plane-parallel H II re- 
gion illuminated by the downstream EUV photon 
field. The density in this region is much higher than 
the pre-ionized region ahead of the shock, so that 
the effective ionization parameter in this zone is 
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Fig. 11. — The oxygen ionization structure, temperature profile and density profile of models with solar abundance, equipartition magnetic 
field for shock velocities of v s = 200, 500, and 1000 km s _1 . The axes are the same as for Figure[3] See the electronic edition of the Journal 
for a color version of this figure. 



much lower, and the H II region is very much thin- 
ner. Eventually the EUV photons are absorbed, 
and hydrogen finally recombines. However, there is 
a cool, partially ionized region behind the main re- 
combination front. This has a temperature of a few 
thousand degrees, and is ionized by the hard pene- 
trating X-rays which undergo Auger processes and 
these lead to the emission of fast electrons which 
in turn produce secondary ionization cascades. 



4.1. Hydrogen ionization structures 



Figures 0] - [7] show, respectively, the effects of chang- 
ing velocity, magnetic field, pre-shock density and abun- 
dance on the structures of the hydrogen ionization frac- 
tion, hydrogen number density, nrj, and the electron tem- 
perature, T c , in the shock and precursor regions. Note 
that in each diagram the time scale changes due to the 
effects of the changing parameters. In each diagram the 
parameters of the model (abundance set, velocity, den- 
sity, and magnetic field respectively) are listed in the top 
right corner for clarity. 

The dominant parameter of the shock (and precur- 
sor) structure and emission is the velocity, as is seen in 
Figure 3J This figure shows our fiducial n= 1.0 cm~ 3 , 
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Fig. 12. — The sulfur ionization structure, temperature profile and c 
field for shock velocities of v s = 200, 500, and 1000 km s . The axes 
for a color version of this figure. 

B/y/n — 3.23 fj,G cm 3 / 2 , solar abundance model with 
two velocities; v s — 200 and 1000 kins" 1 (three includ- 
ing the u s = 500 kms -1 shown in Figure [3]). There is a 
clear increase of both the extent and temperature of the 
radiative zone as the velocity increases due to the greater 
shock strength, leading to the changes in the ionizing ra- 
diation field shown in Figures Q] and [2l This increase of 
both the hardness and luminosity of the radiation field 
leads to the greater extent and ionization in the precur- 
sor region, visible through the increased temperature in 
this figure and the higher ionization species seen in Fig- 
ures I5HT21 The increased radiation field also affects the 
ionization of the post-shock recombination region. This 



lensity profile of models with solar abundance, equipartition magnetic 
are the same as for Figure[3] See the electronic edition of the Journal 

is mollified somewhat however by the density increase as- 
sociated with the faster shock, which leads to a smaller 
increase in the ionization parameter, though the effects 
of the harder radiation are still visible. 

Figure [5] shows models with n=1.0 cm~ 3 , v s = 500 
kms , and solar abundance for differing magnetic fields. 
The top panel shows a model with the magnetic param- 
eter set to B/^n = 0.0001/zG cm 3 / 2 , and the bottom 
panel has B/y/n = 10.0/iG cm 3 / 2 , which can be com- 
pared with the equipartition model with B/ '-y/n = 3.23/iG 
cm 3 / 2 model in Figure [3] In all three cases the precursor 
structure is almost unchanged, as the radiation field and 



14 



Allen et al. 



10 



10' 



22 



10 



20 



c 

0) 

Q 
c 
E 

o 1 
o 



10 1 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 

z 


i | i i i | i 


\ fp^ sir 


- 




- 


' / /'/ 


- 


- DS96 / / 




_ Bn~"=3.23 / Bn"°' 5 =0.0001 

~ i i ■ I i i i I i i ■ I i i 


1 , , , 1 , ~ 



200 



400 600 800 

SHOCK VELOCITY (km s" 1 ) 



1000 



Fig. 13. — Integrated hydrogen column densities for the shock and the precursor structures. The column densities of the shock structure 
are shown for magnetic parameters B/v / n=0.0001, 3.23 and 10. The curve labelled 'PRECURSOR' shows the column density of the 
precursor gas which does not depend on the magnetic parameter. The dashed curves show the shock and precursor column densities for the 
DS96 models, and the dotted curve shows the extrapolation of the column density scaling relation for the shock column density of DS96. 



TABLE 2 
Shock Ionizing Properties 





Luv 


Q(H I) 


Q(He I) 


Q(H II) 


Vion 


T e 


(km s —1 ) 


(ergs cm~ 2 s~ 1 ) 


(km s 1 ) 


(km s — 1 ) 


(km s" 1 ) 


(km s _1 ) 


(K) 


200 


8.041E-4 


209.1 


43.0 


4.9 


214.2 


7033 


300 


3.519E-3 


430.8 


312.1 


83.1 


688.3 


10288 


400 


9.106E-3 


664.3 


598.7 


513.2 


1480.0 


12145 


500 


1.751E-2 


985.2 


1195.0 


602.2 


2318.3 


15435 


600 


2.968E-2 


1158.2 


1628.2 


800.7 


2989.2 


18486 


700 


4.631E-2 


1559.2 


1986.0 


1094.0 


3865.8 


20516 


800 


6.815E-2 


1958.8 


2362.3 


1542.7 


4885.8 


21934 


900 


9.645E-2 


2375.0 


2920.4 


2206.6 


6251.7 


23389 


1000 


1.318E-1 


2740.6 


3262.7 


2786.8 


7325.0 


24760 



luminosity arising from the shock are little affected by 
changes in B (except at extremely high values), as dis- 
cussed in sections $3] and £14.41 The insensitivity of the 
radiation field to B can be seen in this figure through 
the similarity of the radiative zone between the models. 
This similarity, a result of the insensitivity of the temper- 
ature jump to the magnetic field, means that the ionizing 
bremsstrahlung radiation that dominates the energetics 
is practically the same. However, what is sensitive to the 
transverse magnetic field is the density in the recombina- 
tion zone behind the shock front. As discussed in DS96, 
the maximum compression in this post-shock region is 
driven by the ratio of the post- to pre-shock magnetic 
fields. Thus a higher initial magnetic field results in a 
lower density post-shock gas, and therefore, given the 
similarity in radiation field, a higher effective post-shock 
ionization parameter. This results in a significantly dif- 



ferent ionization and emission structure in the post-shock 
gas, leading to the strong diagnostics seen in DS95 and 
in this work 

Figure [5] shows the effect of changing the pre-shock 
density. We plot our fiducial solar abundance models, 
with v s = 500 kms" 1 and B/^n = 3.23 fiG cm 3 / 2 , for 
each of the densities n=0.1, 10, and 1000 cm~ 3 , from top 
to bottom (with the n=1.0 cm~ 3 case shown in Figure 
[3]). Note that the magnetic field, B, is also varied in order 
to keep the magnetic parameter constant for the models 
shown. This figure demonstrates that at a given shock 
velocity, the region influenced by the shock scales with 
density (to the extent that collision de-excitation effects 
are unimportant), proving explicitly that to first order 
the product not is an invariant with changing density, as 
asserted by DS96. 

Changes in the atomic abundances directly affect the 
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Fig. 14. — Total radiative fluxes of shocks as a function of shock 
velocity and various values of the pre-shock densities. The bold line 
is for the standard model having a solar abundance set, a density 
n= 1 cm -3 , and an equipartition pre-shock magnetic field. The 
thinner lines near the bold line show the luminosities of the n= 1 
cm -3 , equipartition models with abundances of twice solar, LMC 
and SMC abundance sets. 

cooling rate. Figure [7] shows the equipartition case, 
= 3.23 ^G cm 3 / 2 , with n=1.0 cm~ 3 and v s = 500 
kms -1 for models using the 2 x Solar, dopita2005, and 
SMC abundances. All three models differ in structure, 
with the expected clear progression in cooling efficiency 
as we increase in metallicity from the SMC abundances 
to the twice Solar. This difference in efficiency leads to 
the changes visible in the spatial and temporal extent 
of the models, in both the precursor and shock regions. 
The changes in metallicity also affect the ionizing radi- 
ation from the shock due to the weakening of the heavy 
element free-bound features. However, as discussed in 
section [3l the effect is minimal. 

4.2. Ionization Structures of Abundant Elements 

The shock model library includes the full fractional 
ionization structure of all elements listed in Table [1] In 
Figures I5HT21 we show the ionization structures of he- 
lium, carbon, nitrogen, oxygen and sulfur for three veloc- 




FlG. 15. — As for Figure \l4\ but showing the precursor total 
radiative fluxes. 

ities (200, 500 and 1000 km s -1 ) of the solar abundance, 
n= 1.0 cm~ 3 models with the equipartition magnetic pa- 
rameter. The associated hydrogen structures are shown 
in Figures [3] and 3J In each figure, the ionization frac- 
tion, X; on , is shown at the bottom (with values < 1), with 
each ionization state labelled and the odd and even ion- 
ization states shown as the solid and dot-dashed curves 
respectively for clarity. Also included in each diagram 
is the hydrogen density, nu, and electron temperature 
structure, T e , to trace the different shock and precursor 
zones. 

Together, these figures elucidate further what has been 
discussed in the previous sections; faster shocks lead to a 
higher temperature, and therefore higher ionization, ra- 
diative zone. This in turn leads to a more luminous and 
harder ionizing spectrum as can be traced within the pre- 
cursor zone. Similarly this figure allows one to trace the 
more complex ionization structure of the post-shock re- 
combination region, affected both by the increased ion- 
izing spectrum and increased density that results from 
faster shocks. 

However the true strength of these figures, and the 
library in general, is their diagnostic ability, allowing the 
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Fig. 16. — Line strength relative to H/3 for 16 strong emission lines between the Lyman Limit and 26fim as a function of shock velocity. 
These are computed for the shock component only, and are for the fiducial model with solar abundance, n= 1.0 cm -3 and B= 3.23 /J.G. 
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reader to trace exactly where each ionization state arises, 
which is especially interesting given the X-ray ionized 
regions that occur far from the ionization front in both 
the precursor and post-shock regions. 

4.3. Column Densities 

The column densities of the shock and precursor com- 
ponents of the models inevitably increase with shock ve- 
locity, and can be as large as iV(H)=4xl0 21 cm~ 2 in the 
shock itself for the highest velocity solar models. The 
precursor column densities range from between 3 and 20 
times larger than the corresponding shock component, 
with this ratio decreasing for the higher velocity shocks. 

Figure [13] displays the integrated H II column den- 
sity of the shock as a function of shock velocity. This is 
shown for the solar abundance n= 1.0 cm -3 models with 
equipartition magnetic parameter, as well as the for max- 
imum and minimum magnetic parameters of 0.0001 and 
10.0. The dashed line shows the scaling relation given 
by DS96, and the dotted line shows the extrapolation of 
this beyond shock velocities of 500 km s _1 . The diagram 
shows a weak dependence of the column density on the 
magnetic parameter, and a decrease in the slope beyond 
shock velocities of 700 km s _1 . 

The column density of the precursor does not depend 
on the magnetic parameter. All solar abundance, n=1.0 
cm~ 3 models fall on the same curve, marked 'PRECUR- 
SOR' in Figure [T3l There is however a systematic dif- 
ference in the precursor column density compared to the 
DS96 models, with the newer models having a ~10% 
decrease in column density, consistent with the slightly 
lower fraction of ionizing flux compared with total flux 
produced in the cooling plasma in the newer models. 

Tables of the column densities for all species for every 
library model are available on-line (see section [5]). These 
columns are a function of both the hydrogen column den- 
sity and the ionization structures seen in the previous 
figures. Example tables giving the column densities of 
the shock and precursor components of the n= 1.0 cm~ 3 
solar abundance models at velocities of 200, 500 and 1000 
km s _1 models are presented in Tables |4]and 5. 

4.4. Radiative fluxes of Shock and Precursor 
Components 

A fully-radiative shock will, by definition, radiatively 
dissipate all of the energy flux through the shock. Thus 
the total radiative flux, Lxot, is equal to 0.5pu 3 . DS96 
provided scaling relations for the total luminosity of a 
shock, and showed a similar scaling exists for the total 
luminosity of the precursor emission. This follows be- 
cause, in the precursor, all of the radiation emitted in 
the pre-shock direction is eventually processed. 

In Figures RBI and [T51 we show the total radiative fluxes 
generated in the shock and in the precursor respectively, 
in the MAPPINGS in models with solar abundance and 
equipartition magnetic fields. Here we cover the full 
range of pre-shock density (0.01 — 1000 cm -3 ). To il- 
lustrate the relatively minor adjustments to the total ra- 
diated flux due to different abundances we also show the 
luminosities for the models with 2x Solar, Dopita2005, 
SMC and LMC abundance sets for n= 1.0 cm~ 3 . Chang- 
ing the magnetic field from the equipartition value pro- 
duces too small an effect on the total radiative flux to be 
visible on the scale of these plots. 



5. SPECTRAL SIGNATURES OF THE SHOCK AND 
PRECURSOR GAS 

While the total luminosities of the shock effectively 
depends only upon the density and the shock velocity, 
the detailed emission line spectra depend strongly on 
the physical and ionization structure of the shock. This 
is determined primarily by the shock velocity and the 
magnetic parameter. However, at the higher pre-shock 
densities (and at lower values of the magnetic parameter) 
the density close to the photoionized tail and the recom- 
bination zone of the shock becomes sufficiently high for 
collisional de-excitation of forbidden lines to become im- 
portant. In these circumstances, the emission line spec- 
trum of the shock becomes dependent upon the density 
as well. 

Shocks are characterized by regions of high electron 
temperature and ionization state. As a consequence, 
they display a rich spectrum of collisionally excited UV 
lines. The shock velocity controls the shape of the ioniz- 
ing spectrum produced by the shock, and the magnetic 
parameter controls the effective ionization parameter in 
the photoionized tail of the shock. Higher ionization pa- 
rameters yield a higher mean ionization in the plasma, 
so therefore faster shocks give an (optical) spectrum of 
somewhat higher ionization. 

The emission line ratios (with respect to H/?=l) 
of a selection of lines for the shock, precursor and 
shock -{-precursor components of the solar abundance, 
n=1.0 cm -3 models are shown separately in Figures \Wl 
IT51 as a function of shock velocity. 

As discussed before, the precursor emission can be con- 
sidered as an H II region with high ionization parameter. 
In this region strong cooling lines, like [O III]A5007, are 
strong and generally increase in strength with shock ve- 
locity as the spectrum produced by the cooling plasma 
becomes harder, and its flux increases. Eventually, as the 
shock velocity is increased, the average state of ioniza- 
tion in the precursor region becomes high and the relative 
intensity of the lines produced by lower ionization turn 
over, flattening and even decreasing at the higher veloc- 
ities. For example, in Figure IT71 [O III]A5007 can be seen 
to flatten out around 600 km s _1 , and even decrease at 
1000 km s" 1 . 

In addition the emission line ratios of 40 strong lines 
are tabulated for shock, precursor and shock+precursor 
in Tables 6, 7, and 8 respectively. Complete tables of 
emission line ratios of all lines, for all models are available 
electronically ( see section [5]) . 

5.1. Line Ratio Diagrams 

Diagnostic emission line ratio diagrams are well estab- 
lished as a powerful and practical way to investigate the 
physics of emission line gas. At a basic level they can 
be used empirically to identify the key excitation mecha- 
nism, and divide active galaxies from their star-formation 
dominated counterparts. At a deeper level, grids of the- 
oretical models can be used to determine chemical abun- 
dances and derive physical parameters such as electron 
temperature, density, ionization parameter or shock ve- 
locity. Various sets o f diagnostic diagrams have been 
proposed and used - Baldwin et alj (|1981) (optical), 
iVeilleux fc Osterbrockl (119871) (on tical-IR) [ Alien et all 
(j!998h Cnear-UV) and iGroves et all (|2004bh (UV, opti- 
cal, IR). These 2-D diagrams remain useful alongside 
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As Figure \TE\ but for the precursor component of the models only. 
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multi-dimensional methods for comparing data to mod- 
els because the ratios can be chosen to be independent of 
reddening, and tuned to be sensitive to pa rticular phys- 
ical p roperties of interest. For example, iDopita et al.l 
(2006) have provided theoretical strong emission line di- 
agnostics for ensembles of H II regions to enable more 
accurate measurement of che mical abundances in unre- 
solved starburst galaxies, and lKewlev et a l. (2006) com- 
bined standard line ratio diagrams with velocity disper- 
sion measurement to help refine the classification of ac- 
tive galaxies. 

As part of the shock model library available online, we 
include interactive plotting tools to plot diagnostic dia- 
grams, such as the examples presented in the following 
sections. These tools enable a much clearer picture of 
these diagnostics than possible within the paper, espe- 
cially with those grids that twist and fold upon them- 
selves in the 2-D diagrams. In addition, many other di- 
agnostic ratios than those presented here can be formed 
to gain insights into what diagnostics can be used when 
only a limited sample of lines are available. 

5.2. Optical Diagnostics 

As a first example, we plot the new shock model grids 
on a number of the standard optical diagnostic dia- 
grams to reveal the general shape of the model grids, 
and the sensitivity of the diagnostics to the parame- 
ters of shock velocity, magnetic parameter and chem- 
ical abundance set. Figure [19] shows the perhaps the 
most frequently used of the diagnostic diagrams, the 
IVeilleux k Osterbrockl (fl98l plot of [O III]A5007/H/3 
versus [N IIlA6583/Ha , also comm only known as the 
BPT diagram (|Baldwin et al.lll981f) . For clarity, the 
shock and shock+precursor grids for the solar abundance 
n= 1 cm~ 3 models are plotted separately but on the 
same scale, as is the case in most of the following dia- 
grams, due to the overlap of the two grids. These two 
grids represent the physical extremes of having little or 
no gas ahead of the shock (shock-only), and having an 
extensive, radiation bounded, precursor region ahead of 
the shock. As in DS95, the two grids show the range 
in shock velocity and magnetic parameter, as labelled in 
the diagram. The lines of constant velocity are plotted 
at 50 km s _1 intervals as the thin lines. The lines of con- 
stant magnetic parameter are plotted as the thick lines 
in the grid, with the value of the magnetic parameter as 
marked. For a clearer picture of this and all other diag- 
nostic diagrams in the paper we provide color figures in 
the online journal. 

This figure reveals a similar picture to that in DS95 
(their Fig. 2b), with the shock-only and shock+precursor 
grids occupying different regions of the diagram, a result 
of the different processes discussed in the previous sec- 
tions. Both grids shows broad steps in velocity up to 
a velocity of ~ 550 km s -1 , at which the grids begin to 
compress and turn upon themselves, becoming degen- 
erate. However, for the shock-only grid an increase in 
velocity causes an increase in the [N II] /Ha ratio, while 
for the shock+precursor models in results in an increase 
of the [O III]/H/3 ratio. The situation is reversed for an 
increase of the magnetic parameter, causing the spread 
observed in the grids, though the shock+precursor mod- 
els tend to be degenerate in this parameter at high values. 
The spread of the grids and their separation indicate the 



diagnostic power of these figures, allowing the determi- 
nation of the physical parameters in known shock-ionized 
plasmas. 

Improvements to the MAPPINGS in code compared 
to the MAPPINGS II code used by DS96 means that 
there are some significant and systematic differences 
compared to the DS95 grids. Figure [20] shows a di- 
rect comparison of the DS95 shock and shock+precursor 
model grids overlaid on our new solar abundance shock 
and shock+precursor models for the same diagnostic di- 
agram. The improvement of our new models in terms of 
parameter space is clear, with DS95,96 models limited to 
500 km s _1 in velocity space, and only covering magnetic 
parameters of 0-4 fj,G cm 3 / 2 . The shock component of 
the models agree reasonably well over the common pa- 
rameter ranges, although appears somewhat offset. How- 
ever the shock+precursor grid shows a difference of up 
to 0.5 dex in the [O III]A5007/H/3 and [N II]A6583/Ha 
ratios. We note again that on this diagram the new 
shock+precursor models turn over at high shock veloc- 
ity. One consequence of this is that some active galax- 
ies such as high-redshift radio galaxies, which have been 
confirmed as being excited by shocks and which were in- 
terpreted as matching 400-500 km s _1 shock+precursor 
models, may in fact be characterized by higher shock ve- 
locities. This would help resolve the mis-match between 
radial velocity dispersions and shock vel ocities inferred 
by lin e ratios as inferred, for example, bv lReuland et al] 
(|200l . 

Figure [21] shows the same diagram for the different 
chemical abundance sets which we have used. Each grid 
covers the same range in shock velocity and magnetic pa- 
rameter. In both shock and shock+precursor models the 
[O III]A5007/H/3 ratio is only weakly affected by abun- 
dance, as the increased temperatures within the models 
from the reduced cooling counteracts somewhat the ef- 
fects of the drop in oxygen abundance on this ratio (see 
Figure [7]). The [N II]A6583/Ha ratio however changes 
significantly with abundance, predominantly as a conse- 
quence of the larger relative abundance differences in ni- 
trogen due to its secondary nucleosynthesis component. 
In the shock models, the [O III]/H/3 does not change 
significantly in range either, but the [N II] /Ha ratio 
becomes compressed, meaning the diagrams are totally 
degenerate in velocity at low metallicity, predominantly 
due to the changes photoionized and recombination re- 
gions of the shock with metallicity. The case is sim- 
ilar for the shock+precursor, with the grids becoming 
generally more "compressed" with decreasing metallic- 
ity, and decreasing the diagnostic power of the diagrams. 
Note that for each abundance set the shock-only and 
shock+precursor models still tend to be discrete grids, 
and do not significantly overlap. 

In Figure [22] we investigate the effect of changing 
the p re-shock densities on this IVeilleux k Osterbrockl 
l|l987h diagnostic. For the solar abundance grid, we dis- 
play the results for the two extreme values of density 
n= 0.01 cm~ 3 and n= 1000 cm' 3 for both the shock 
and shock+precursor models. Note the change in scale 
between this figure and Figure \T§\ necessary due to the 
changes in the grids. 

As discussed in section 12.11 the models with varying 
density were computed not only with the standard mag- 
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Fig. 18. — 



As Figure IT6l but for the the full shock+precursor model. 
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Fig. 19. — The classic Vcillcux & Ostcrbrock (1987) [O III] A5007/H/3 versus [N II]A6583/H« diagnostic diagram for the solar abundance 
n= 1 cm -3 models. The left figure shows the shock-only models, and the right figure shows the shock+precursor models. The grid is 
comprised of lines of constant magnetic parameter shown with a bold linestyle, and lines of constant shock velocity shown as the thin lines. 
The lines of constant shock velocity are shown at 50 km s _1 intervals (and are color coded with increasing red intensity for higher shock 
velocities in the electronic edition). The velocity sequences are shown for all values of the magnetic parameter (and are color coded with 
increasing blue intensity for higher magnetic parameters in the electronic edition). Both the shock-only and shock+precursor grids are 
shown for the shock velocities of 200-1000 km s —1 , for which all the models have a fully ionized precursor. See the electronic edition of the 
Journal for a color version of this figure. 

models in Figure [2"2l for completeness, meaning that the 
grids sample more finely and extend further both the low 
and high values of the magnetic parameter than seen in 
Figures HUT] 

The low-density model grids are in most respects ex- 
actly the same as the n= 1 cm -3 grids seen in Figure 
[T9l with the only difference arising due to the extended 
B/y^ range sampled. This similarity is not surprising, 
as while the luminosity decreases linearly with the de- 
creased density, as seen in Equation [5] the ionization 
parameter is inversely proportional to density and there- 
fore remains approximately the same in both the shock 
recombination region and the precursor region. One in- 
teresting thing to note is how degenerate the grids be- 
come upon the introduction of the very high magnetic 
parameters. 

However, at high pre-shock densities there are clear 
changes in both the shock and shock+precursor grids. 
In the shock-only grid the low magnetic parameter 
models extend to lower values in the both [O IIl]/H/3 
and [N II] /Ha ratios by up to 1 dex, while in the 
shock+precursor grid only the [N II] /Ha ratio is reduced. 
This is caused by collisional de-excitation of the forbid- 
den lines in the photoionized and recombination regions 
of the shock, where the electron densities are very high. 
The shock+precursor does not change significantly in 
[O III]/H/3 as most of the emission of the [O III] line 
arises from the precursor, which is not high enough den- 
sity to collisionally de-excite this transition. There is less 
of a difference between the grids at high magnetic pa- 
rameters and at low to intermediate velocities, because 
lower velocities lead to smaller compression factors, and 
because magnetic pressure support helps to lessen the 
degree of gas compression in the tail end of the shock, 
meaning that the high density needed to collisionally de- 
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Fig. 20. — Comparison of the new MAPPINGS III models with 
the DS96 on the [O III] A5007/H/3 versus [N II]A6583/Ha diagram. 
The DS96 models are shown in black, they are labelled with the 
shock velocity in km a -1 , and include velocity sequences for mag- 
netic parameters of 0, 1, 2 and 4 /xG cm 3 / 2 which increase from left 
to right in both the shock-only (lower) and shock+precursor (up- 
per) grids. The new MAPPINGS III models are shown in greyscale 
and cover the same range in shock velocity and magnetic field as 
shown in Figure [19] Note that the new diagrams show systematic 
changes with respect to the earlier ones, and that the new mod- 
els show that these ratios become insensitive to shock velocity for 
v g > 500 km s _1 . 

netic parameter set, but also a range of magnetic fields to 
allow for the comparison of both B and B/ym" between 
the model grids. The range of both B and B/ym" for 
these models are given in Table [3] We include all these 
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Fig. 21. — Shock and shock+precursor model grids for the five different atomic abundance sets used in the shock model library for the 
[O III]/H/3 versus [N II] /Ha diagram. The shock-only models are shown in (a), and the shock+precursor models are shown in (b). All 
model grids are shown for d s 200—1000 km s — 1 , with lines of constant shock velocity drawn at intervals of 50 km s — 1 . Each grid is labelled 
with the abundance set used, and generally moves from left to right as the total metallicity increases. See the electronic edition of the 
Journal for a color version of this figure. 

sion in the photoionized tail of the shock leads to a higher 
ionization parameter in this zone, and greater dominance 
of [O III] in the photoionized zone, which has an elec- 
tron temperature close to 10000K. As a consequence, the 
electron temperature decreases and the [O III]A5007/H/3 
ratio increases with increasing shock velocity. Together 
this leads to the fan shape observed in the diagram. The 
shock+precursor grid is dominated by the bright [O III] 
emission of the precursor region, and as the precursor is 
not sensitive to the magnetic parameter, this grid covers 
a reduced area of the parameter space in the diagram. In 
terms of velocity, the curve increases in [O III]A5007/H/3 
and decreases in i?om up to a velocity of ~ 500 km s -1 , 
then turns over and begins to increase in temperature 
again. The shock+precursor grid at velocities less than 
^450 kms -1 and the high magnetic parameter shock- 
only grid overlap, meaning this diagnostic is degenerate 
in this range. This diagram is very similar to the DS95, 
with only the extension to higher velocities and mag- 
netic parameter range differentiating it. However, this 
similarity also indicates that the well-known "Tempera- 
ture Problem" , discussed in DS95 and many other pa- 
pers, still exists. A comparison of our models with the 
observations used in DS95 or larger datasets like SDSS 
show that the i?oiu ratio is too low, and even different 
abundance sets are unable to fully solve this issue. 

5.3. UV Diagnostics 

In I Allen et al.l (|1998f ) we emphasized the use of UV line 
ratio diagrams for the discrimination between shocks and 
photoionization models for the NLR. These diagrams use 
relatively bright emission lines and reddening-insensitive 
ratios. One of the most useful diagrams involves the 
various ionization stages of carbon C II]A2326/ C III] 
A1909 vs. C IV A1550 / C III] A1909. This diagram is 
shown in Figure [55] with the new solar abundance n= 1 
cm~ 3 shock models. There is significantly more overlap 



excite these transitions is not rea ched. 

Ano ther commonly used IVeilleux fc Osterbrockl 
(1987) diagnostic is the [O III]A5007/H/3 versus 
[S II]AA6716,6731/Ha diagram. In Figure [23] we plot 
the new grids for this diagram, with the shock and 
shock+precursor models plotted separately for clarity, 
as the overlap between the grids makes distinguishing 
them difficult. As with the previous diagnostic diagram, 
both the shock and shock+precursor models turn 
over at high velocities, becoming twisted and with 
little separation between the velocities. Similarly the 
magnetic parameters show a similar range in spread in 
both models, and become more spread with increasing 
velocity in the shock-only model. The shape of the 
shock+precursor grid on this diagram is again rather 
different compared to DS95. Like in Figure [20l the 
new shock+precursor model covers a wider range in 
[O III]/H/3, with a lower value at v s = 200 kms -1 , while 
it appears somewhat offset to higher values of [S II] /Ha, 
though not significantly. It also does not appear to show 
the degenerate nature see in DS95 (their Figure 2a). 
The shock models are harder to compare due to the 
wider parameter range in the model, but in general show 
the same features of increased [S II] /Ha with velocity, 
with the turn over around ^500 kms -1 . The greater 
range in magnetic parameter of the new models also 
means that the spread in [O III]/H/3 is much greater at 
higher velocities as well. 

Figure [24] is a plot of the well-known temperature- 
sensitive ratio [O III]A4363/[0 III]A5007 (also known 
as i?om) versus the [O III]A5007/H/3 ratio, with again 
the shock-only and shock+precursor spectra plotted sep- 
arately. Considering first the shock-only grid, at low 
magnetic parameter the electron temperature in the 2+ 
zone is consistently high, with only the [O III]A5007/H/3 
ratio decreasing with increasing shock velocity. However, 
as the magnetic parameter rises, the decreased compres- 
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Fig. 22. — Shock and shock+precursor model grids of solar abundance models for two different pre-shock densities on the [O III] A5007/H/3 
versus [N II]A6583/Ha diagram. The top figures show the high density case and the lower figures the low density case, with a) The n= 1000 
cm -3 shock model grid, b) the n= 1000 cm -3 shock+precursor models, c) The n= 0.01 cm -3 shock model grid, and d) the n= 0.01 cm -3 
shock+precursor models. The grids show shock velocities over the range d s =200-1000 km s -1 with lines of constant magnetic parameter 
shown as thick lines (and colored with increasing blue intensity for higher magnetic parameter in the electronic edition). The thin lines 
represent constant shock velocity and these are plotted for the full range of magnetic parameters for each density as listed in Table 2 (and 
are shown colored with increasing red intensity for higher shock velocities in the electronic edition). See the electronic edition of the Journal 
for a color version of this figure. 



between the new shock and shock+precursor grids than 
with the DS96 models because the higher magnetic pa- 
rameter shock models form a fan-like grid on these axes, 
which almost completely encompasses the region covered 
by the shock+precursor models. The shock+precursor 
grid is twisted on these axes making it multi-valued at 
most positions, hence not good for assessing shock pa- 
rameters. This diagram does however remain of great 
use for separating shock and photoionization models. 

Another useful diagram combines the UV carbon ra- 
tion of C III] A1909 / C II] A2326 wit h [Ne III] A3869 / 
[Ne V] A3426. This ha s been used by iBest et all (|2000l ) 
and llnskip et a l. (2002) to identify shocks and photoion- 
ization in the emission line gas of 3CR and 6C radio 



galaxies. They show that the ionization state of the 
gas varies with radio size such that large radio sources 
(> 120 kpc) are consistent with AGN photoionization, 
while smaller sources are consistent with shocks associ- 
ated with the expanding radio source. Furthermore the 
extreme gas kinematics in the smaller radio sources, and 
detailed con sideration of t he en ergetics and observability 
of shocks in llnskip et alj (|2002D supports the interpreta- 
tion of shock excit ation. 

As described in llnskip et alj (|2002f ). there is a con- 
tinuous sequence of objects which fall in between the 
shock dominated and photoionization dominated groups. 
In these, both mechanisms are likely to play a role and 
the balance between shocks and photoionization is may 
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Fig. 23.— The lVeilleux fc Ost crbrock (1987) plot of [O III]A5007/H/3 versus [S II]AA6716,31/Ha diagram for a) the shock and b) the 
shock + precursor models with n= 1 c m" 3 and solar abundance. The range and step size of the shock velocity and magnetic parameter 
for these grids are the same as in Figure [19] See the electronic edition of the Journal for a color version of this figure. 
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Fig. 24.— The temperature sensitive [O III]A4363/ [O III]A5007 ratio (-Roill) versus the [O III] A5007/H/3 line ratio diagram for the 
solar abundance n= 1 cm -3 models. Panel (a) shows the shock grid, and panel (b) shows the shock+precursor grid. The range and step 
size of the shock velocity and magnetic parameter for these grids are the same as in Figure [19] See the electronic edition of the Journal for 
a color version of this figure. 



be linked to the radio source size. As the definition of 
the shock and photoionization groups relies on the cov- 
erage of the respective model grids, it is important that 
these grids are complete. In Figure[26]we reconstruct the 
same C III] A1909 / C II] A2326 versus [Ne III] A3 869/ 
[Ne V] A3426 diagr am as used bv lBest et ail (|2000T ) and 
llnskip et all {2002) using the data tabulated in those pa- 
pers. We overlay the new shock model grids to demon- 
strate how these relate to, and possible help explain, 
these previously-made conclusions. For comparison, we 
also plot the set o f new and updated AG N photoioniza- 
tion models from iGroves et al.l (l2004bf) as well as the 
A M /j models from iBinette et al. I (119961) that combine 



matter- and ionization-bounded clouds to create the ob- 
served sequence. 

The shock and shock+precursor models are those for 
new solar abundance n= 1 cm~ 3 . As in the previous 
diagram, the shock models define a fan shaped grid, 
and the shock+precursor grid turns over in both ra- 
tios. With the wider range in shock velocity and mag- 
netic parameter the new models cover a great er region of 
this d iagra m than the DS96 m odels used by iBest et ahl 
(|2000h and llnskip et~aT] (|2002D . Indeed a number of 6C 
sources which have C III] / A1909 C II]A2326 ratios in- 
termediate between the main shock and photoionization 
groups, fall within the higher velocity shock+precursor 
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Fig. 25. — The ultraviolet carbon line ratio diagnostic plot for the solar abundance n= 1 cm -3 shock (a) and shock+precursor models 
(b). The shock only grid enables both the shock velocity and the magnetic parameter to be independently determined, but the shock + 
precursor grid is mult i-va lued at many positions. The range and step size of the shock velocity and magnetic parameter for these grids are 
the same as in Figure 1191 See the electronic edition of the Journal for a color version of this figure. 
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Fig. 26.— The C III]A1909/C IIJA2326 with Ne IIIA3869/Ne V]A3426 diagram. The shock model grid is shown in light gray (cyan 
in the electronic edition) and is labelled with the shock velocity. The shock+precursor grid is shown in dark grey (green in the electronic 
edition ), and displays a twisted shape on these axes, and also overlaps the shock grid. The dusty and dust-free models of lGroves et al.l 
Il2004bl ) are also shown with the ionization parameter increasing to the left of the diagram before the turn-over in these curves. The A M yj 
models from Bincttc ct al. ( 1996) are shown a s a do tted line with the + symbols indicating increments of 0.25 dex. The data points are 
those from Best ct al. (2000) and Inskip ct al. (2002) as indicated in the text. 
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Fig. 27. — IR diagnostic diagram of [Ne V]14.3/im/[Ne II]12.8^m versus [Ne III]15.5/^m/[Ne II]12.8/im. Panels (a) and (b) show 
respectively the shock and shock+precursor model grids for the solar abun dan ce n= 1 cm" 3 models. The range and step size of the shock 
velocity and magnetic parameter for these grids are the same as in Figure [19] See the electronic edition of the Journal for a color version 
of this figure. 
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Fig. 28.— IR diagnostic diagram of [Ne V]14.3/mi/[Ne II]12.8^m versus [O IV]25.9/tm/[Ne II]12.8/im. Panels (a) and (b) show 
respectively the shock and shock+precursor model grids for the solar abun dan ce n= 1 cm -3 models. The range and step size of the shock 
velocity and magnetic parameter for these grids are the same as in Figure IT9l See the electronic edition of the Journal for a color version 
of this figure. 

grid (6C 1017, 6C 1256, 6C 0943, 6C 1019). 

The dust-free and dusty photoionization models from 
iGroves et al.1 (|2004bf ) demonstrate the other extreme of 
ionization, where the emission arises totally from gas ex- 
cited by the ionizing radiation emitted by the accretion of 
gas onto the central black hole of the AGN. The dust-free 
models represent MAPPINGS in photoionization models 
that have been calculated for a sequence of ionization pa- 
rameters and densities. The 1Z Q , a=-lA, nH=10 2 -10 4 
models shown here behave in a similar way t o the MAP- 
PINGS II models described in lAllen et alj (1998) (and 



used bv lBest et ail (|2000D and llnskip etTaTl (|2002D V The 
dusty models incorporate the effect of radiation pressure 
and result in a stagnation of the ionization parameter 
at high values, offering an explanation for the similar- 
ity of Seyfert NLR spectra. Figure [26] shows that the 
dusty photoionization models can produce line ratios as 
observed in the larger 3CR and 6C radio sources. The 
dusty photoionization model plotted here is the 1Z Q , a—- 
1.4, nn=10 2 and nH= 10 4 model, and a d iscu ssion of these 
model s in terms of the lBest et al.l (|2000f ) and llnskip et al.l 
( 2002) data, as well as a full grid of such models on the 
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Fig. 29.— IR diagram of [S IV]10.44/im/[S III]18.7Atm 
versus [Ne III]15.5/^m/[Ne II]12.8/im showing the shock and 
shock+precursor grids of the solar abundance n= 1 cm" 3 mod- 
els, as labelled. The range and step size of the shock velocity and 
magnetic parameter for these grids are the same as in Figure [19] 
Sec the electronic edition of the Journal for a color version of this 
figure. 
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Fig. 30.— IR diagram [Fe II]26.0/im/[O IVl25.9^m versu s 
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same diagram can be found in lGroves et"aL1 (|2004bf) . 

5.4. IR Diagnostics 

The mid- and far-IR emission from galaxies is domi- 
nated by the emission of dust that is heated by UV radi- 
ation. IR spectra contain PAH and dust features and are 
also rich in atomic fine structure emission lines. These 
spectral features, as observed by ISO and Spitzer offer 
a wealth of information for studying the nature of the 
circumnuclear dust, and the contributions of the AGN, 



1.5 
1.0 

X 

\ 

1^ 
o 

o 

s 0,0 

f -0.5 
-1.0 



-1.5 



-|— i— I— I— I— I— I— I— |—r- 

SH0CK+PRECURS0R 



200 




-2.0 



-1.5 



-1.0 -0.5 0.0 
log( [Nil] X6583 / Ha ) 



0.5 



1.0 



Fig. 31. — Comparison of the shock and shock+precursor models 
to SDSS line ratios on the [O III] A5007/H/3 versus [N II]A6583/Ha 
diagram. The model grids are those with 2 X solar abundance. The 
dot ted and dashed curved lin es represent the Kcwlcv et al.l (l2001aD 
and IKa uffmann et al. (2003) classification lines respectively which 
divide H II regions (lower left) from the region occupied by AGN. 
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Fig. 32. — Comparison to SDSS line ratios to the 2xsolar 
abundance model gr ids. The thin dotted line represents the 
IKewlev et al.l (l2001aD starburst / AGN classification line, while the 
thick straigh t dotted line repres ents the Seyfert-LINER division 
described by Kcwlcv et al. (2006). 



s hocks and starbursts to the total IR emission. 

iGenzel et al.l (|1998l) introduced infrared diagnostic di- 
agrams using ratios of [Ne V]14.3/im / [Nc II]I2.8/im 
and [O IV]25.9/im / [Ne II]12.8^m to investigate star- 
bursts, ULIRGs and AGNs . They used these line ratios, 
combined with PAH strengths to place limits on the per- 
centage contribution of AGN and starburst contributions 
to the IR emission of U LIRGs. Mid-IR di agnostics for 
LINERs are described bv lSturm et al.l (|2006f ). where they 
use the high ionization [O IV]25.9/im and [Ne V]14.3/im 
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Fig. 33. — Comparison to SDSS line ratios to the 2x so- 
lar abundanc e mode l grids. The thin dotted line represents the 
Kcwlcy ct al. (2001a) starburst/AGN classification line, while the 
thick straigh t dotted line repres ents the Seyfert-LINER division 
described bv lKewlev et al.l H 20061) . 

lines to identify the presence of AGN in ~ 90% of LIN- 
ERs. They confirm the differences in properties between 
IR-faint and IR-luminous LINERs, and identify the need 
to disentangle the various stellar, H II region and AGN 
p rocesses at work in these objects. 

iGroves et al.l (|2006f ) investigated the IR emission of the 
NLR distinct from the emission of the torus, showing 
that the NLR emission can contribute up to ~ 10% of 
the IRAS 25/im flux. They emphasize the fact that high 
ionization lines like [Ne V]14.3/im arises only in the NLR, 
and their diagnostic diagram which utilizes only IR lines 
of neon provides a very useful indicator of AGN and star- 
burst contributions. 

Here we present a set of IR diagnostic diagrams 
drawn from these previous works, overlaid with the 
solar abundance shock and shock + precursor model 
grids. Figure [27] ([Ne V]14.3/im / [Ne II]12.8Atm versus 
[Ne III]15.5/^m / [Ne II]12.8/^m) combines three differ- 
ent ionization stages of neon removing any abundance 
dependence. Figure |2"51 [Ne V]14.3^m / [Ne II]12.8/xm 
versus [O IV]25.9/im / [Ne II]12.8/im, utilizes two high 
ionization species of O IV and Ne V whose ratios to 
Ne II provide strong discriminants of AGN versus star- 
burst processes. Figure [29] displays the model results 
for [S IV]10.44^m / [S III]18.7^m versus [Ne III]15.5Atm 
/ [Ne Il112.8u m. 

iLutz et all (|2003f ) used the ratios of [Fe II]26.0/im / 
[O IV]25.9/Lfm versus [O IV]25.9^m / [Ne II]12.8/im to 
distinguish between gas ionized by early-type stars, AGN 
and shocks. ISturm et alj (]2006[ ) showed how the same di- 
agram separates starburst galaxies, Seyfert galaxies and 
supernova remnants. Figure [301 shows how the shock and 
shock + precursor model grids form relatively tight shock 
velocity sequences in this diagram. The slope of these se- 
q uences is simi l ar to the distribution of observed ratios 
in ISturm et al.l (|2006l ) but the models predict systemati- 
cally higher ratios than observed, except for some LINER 
objects which fall in the region of the lower velocity shock 



model grid. 

5.5. Comparison with SDSS Observations 

We now compare the new shock and shock+precursor 
models to emission line ratios of AGN and star-forming 
galaxies observed in the Sloan Digital Sky Survey 
(SDSS). To do this we have adopt ed the sample of n arrow 
emission line galaxies compiled bv lHao et al. (2005). The 
classification of the sample into broad- and narrow-line 
A GN, and st a r-form ing galaxies is considered in detail 
in lHao et all (|2005l ) . Here we choose to use the com- 
plete narrow emission line galaxies sample of ~ 42000 
sources, and us e the AGN/star-form i ng galaxy separa- 
tion s cheme s of |Kewlev et al.l (l2001al ) , iKauffmann et all 
( 2003), and iKewlev et al.l (|2006l ) in order to emphasize 
where our model grids lie with respect to these different 
classes of objects. 




-2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 
log( [Ol] *6300 / Ha ) 



Fig. 34. — Comparison to SDSS line ratios to an oxygen exci- 
tation diagnostic in which the regions identified as being excited 
by H II regions, Seyferts and LINERs are labelled. Note that 
shock-only spectra may characterize some LINERs, but that the 
Seyfert galaxies are not well fit by shock+precursor models on this 
particular diagram. 

Figu re | 31 | shows the first of the famil- 
iar IVeilleux fe Osterbrockl (|1987l ) diagnostics; 
[O III]A5007/H/3 versus [N II]A6583/Ha. We have 
overlaid as a density plot the line ratios observed in the 
SDSS the narrow emission line sample. 

The distribution of emission line galaxies shows two 
main branches. The star forming galaxy branch sweeps 
in a curve showing relatively small scatter from upper left 
to lower right. This is largely an abundance sequence, 
with abundances increasing towards th e lower right 
(iDopita et all l2000l : IKewlev et all l2001bl : iDopita et all 
2006?). The AGN are mostly distributed in an arm which 
extends from the base of the star-forming sequence up 
towards the upper right of the diagram. AGN gener- 
ally have higher values of [N II]A6583/Ha, and the dis- 
tribution extends to higher values of [O III] A5007/H/3. 
The dotted and dashed curved lines represent the 



The MAPPINGS III Shock Model Library 



29 



iKewlev et all (|2001al) and IKauffmann et al.l (120031) clas- 
sification lines respectively The IKewlev et a" (2001a) 
classification is based on the theoretical maximum line 
ratios possible by pure stellar photoionization. Galaxies 
above this line are most likely dominated by AGN. The 
IKauffmann et all (|2003f ) line is a purely empirical divid- 
ing line between pure star-forming galaxies, and Seyfert- 
HII composite objects. 

For comparison, we show in Figure [3T] the shock and 
shock -(-precursor grids for the 2 x solar metallicity mod- 
els. This cho i ce of meta llicity is driven by t he work of 
iGroves et all (|2004f ) and IGroves et all (|2006| ). who find 
that super-solar metallicity photoionization models best 
reproduce the observed narrow-line ratios in AGN, and 
the work of lKauffmann et all (|2003l ). who find that AGN 
are typically hosted in galaxies more massive than 10 10 
M Q , and therefore likely t o contain high metallicity gas 
(e.g. iTremonti et~aTll2004D 

Figure [31] reveals that both the shock and 
shock -(-precursor models are generally located in 
the AGN region o f this diagram, lying mostly above 
the IKewlev et al.l (|2001al ) classification curve. The 
shock -(-precursor grid overlaps well wit h the strong 
AGN or Seyfert branch, lying above the IKewlev et al.l 
(2006) LINER/Seyfert dividing line, and extends along 
the branch with increasing shock velocity. The high 
shock velocity portion of the grid extends roughly to the 
limit of the observed distribution before folding over on 
itself at the highest velocities. At lower velocities, the 
shock -(-precursor models extend into the "composites" 
region an d even into the s t arfor mation or HII region 
below the IKauffmann et all (|2003l ) curve at the lowest 
velocities and magnetic parameters. 

The lower velocity shock-only models also overlap the 
AGN branch, but are predominantly located in the 
LINER region, simultaneously extending out and spread- 
ing out (with magnetic parameter) at the higher shock 
velocities. 

The emission -line galaxies form sim i lar d istributions 
on the other two lVeilleux fe; Osterbrockl (jT987) diagnostic 
diagrams; [O III] A5007/H/3 versus [S II]AA6716,6731/Ha, 
and the plot of [O III] A5007/ H/3 versus [O I jA630 0/Ifa 
(see Figures [32] and EH} . The IKewlev et all |2001ah di- 
vision between AGN and H II region excitation is shown 
as dotted lines on these figures. 

In these figures, the A GN branches s h ow a b ifurcation, 
that has been used by IKewlev et al.1 (|2006l ) to distin- 
guish between Seyfert and LINER galaxies, and so pro- 
vide a general classification scheme for AGN host galax- 
ies. The Seyfert-LINER dividing line is shown in Fig- 
ures [32] and [33] as a thick dotted line. The shock and 
shock -(-precursor grids overlap the AGN distributions, 
but tend to fall in the region of the Seyfert-LINER di- 
vision. The shock-only models mostly fall in the LINER 
region of these diagnostic plots, but tend to extend to 
extend to higher values of [S II]AA6716,6731/Ha, and 
[O I]A6300/Ha! than is observed. There is also more 
overlap with the shock+precursor models than in seen 
in Figure I3T1 

Figure M shows the [O III]A5007/[O II]AA3726,3729 
versus [O I]A6300/Ha diagnostic diagram over-plotte d 
with the classification scheme from IKewlev et ahl ([2006). 
This diagram provides the cleanest separation between 
the two AGN branches. In photoionized plasmas, the 



[O III]A5007/[O II]AA3726,3729 is sensitive to the specific 
intensity of the radiation field, and the [O I]A6300/Ha 
ratio to the hardness or spectra l index of th e radi - 
ation field. As pointed out by IKewlev et alj (|2006f ) 
this diagram provides a simple method for classifica- 
tion, but is more sensitive to reddening correction of 
the [O III]A5007/[O II]AA3726,3729 ratio. Once again, 
the shock+precursor grids fall predominantly close to the 
Seyfert-LINER dividing line, and so describe neither the 
Seyferts or the LINERs particularly well. The pure shock 
models are a much better fit to the LINER sequence, 
provided that the shock velocities are not too great. The 
highest shock velocities have too strong [O I]A6300/Ha 
ratios. 

In conclusion, Seyfert galaxies are not well described 
by the shock + precursor models. For these, the 
radiation-pres s ure dominated photoioniz ation models 
(jDopita et all 120021 : IGroves et all l2004al lbT) provide a 
much better description of the spectra. The LINERs, 
on the other hand, fit better to shock only models, and 
it is likely that at least some of these objects are in fact 
shock-excited. A good example of a LINER which is 
known to be shock -excited is the nuclear disk of M87 
(jDopita et al.lll997f) . 

6. THE ON-LINE LIBRARY 

The complete electronic files that comprise 
the MAPPINGS in Library of Fast Shock Mod- 
els are available via the Shock Model Por- 
tal of the MAPPINGS online web pages at 
|http : //cdsweb . u-strasbg . f r/~{}allen/ shock . html[ 
The original MAPPINGS in output files are available 
for each model in the library. We also provide tables of 
emission line ratios, and column densities. These files 
are organized into the various model sets characterized 
by a given chemical abundance and pre-shock density 
and then into velocity sequences for a given magnetic 
field. The files are listed on the online pages using the 
model names shown in Table [3] 

The emission line ratio tables are wavelength ordered 
lists of the flux ratios of all the emission lines calculated 
in MAPPINGS in, and are given with respect to H/?=l. 
Each table contains the ratios for a velocity sequence 
of models, 100-1000 km s _1 , for a given abundance set, 
density and magnetic field. Shock and precursor compo- 
nents are tabulated separately and we also provide ta- 
bles of the emission line ratios for the combination of 
shock+precursor. Each emission line ratio table also in- 
cludes the absolute luminosity of the H/? line in units of 
Logio(erg cm~ 2 s _1 ). 

The column density tables contain the integrated 
model column densities for all of the ionic stages of 
each of the elements in the corresponding abundance list. 
Each table includes the column densities for a velocity 
sequence of models in units of cm~ 2 . Column densities 
for the shock and precursor components are tabulated 
separately. 

In addition to the tabulated model data, we also pro- 
vide programs for accessing the library of models, and 
for generating various plots. These programs are coded 
using IDL 3 and can be used as interactive graphical user 
interface widgets, and also via the IDL command line. 

3 http://www.ittvis.com/idl/ 
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The SHOCKPLOT package allows plotting the models 
on 2-D line ratio diagrams using any linear combina- 
tion of line ratios. SHOCKPLOT enables quick browsing 
through the many grids of models, and can also be used 
as a procedure call from other IDL programs to over-plot 
model grids on observed line ratio data. 

7. SUMMARY 

We have presented an extensive library of radiative 
shock models covering a wide range of shock velocities, 
magnetic fields, densities and abundances. The shock 
model predictions for the ionizing radiation, tempera- 
tures, and luminosities generated by shocks will be ap- 
plicable in a wide range of astrophysical situations. The 
solar abundance models supersede the models of DS96, 
with model code improvements leading to some differ- 
ences in the predictions of nitrogen species, and also 
for models with very low magnetic parameter. Differ- 
ences in the input abundances result in the most sig- 
nificant variations between the grids of models, and the 
range and sampling of the shock velocity and magnetic 
field provides much a more detailed and complete set 
of model predictions than previously available. The ex- 
tension to higher shock velocities is important for the 
analysis of emission line regions of active galaxies where 
shocks of this speed are expected in jet-cloud interactions 
(jSaxton et alJ 12005). At these shock velocities some of 
the commonly used line ratios such as [O III]A5007/H/3 
turn over and do not follow simple extrapolations from 
the DS96 model grids. 

Included as a part of this library, the physical struc- 
ture of the shocks and their precursors give insight into 
how the physical parameters of the shock lead to the re- 
sulting continuum and line emission. We have included 
in this paper examples of these, exploring the range of 
the parameters and demonstrating the effects of these of 
the density, temperature and ionization structure. 

We have presented the model grids on a set of standard 
UV, optical and IR line ratio diagrams, and we have 
compared of the new models to the example data sets of 
radio galaxies and the SDSS sample of narrow emission 
line galaxies. The updated version of the C III]A1909 
/ [C II1 A2326 vs. [Ne III1A 3869 / [Ne V1A34 2 6 diag ram 
used by iBest et al.1 ([20001 ) and llnskip et all ([20021 ) for 



the analysis of emission line regions of radio galaxies, 
shows that some of the 6C sources fall within the new 
higher velocity shock+precursor grid. This supports the 
interpretation that the emission line regions in smaller 
radio sources are excited by shocks. The new library 
of models presented here, and a more complete set of 
emission line observations of these sources will allow a 
much more detailed analysis of the contributions of shock 
and photoionization in these sources. 

Large samples of emission line ratios such as now avail- 
able from the SDSS provide an extremely valuable re- 
source for classifying and analyzing the emission line ex- 
citation mechanisms in active and star- forming galaxies. 
The comparison of the narrow em ission line galaxies sam- 
ple (compiled by lHao et al.l (|2005h ) to the shock models 
shows that shocks do predict line ratios in the observed 
range, and that shocks may indeed provide the best 
explanation for LINER spectra. This library of shock 
models, combined with detailed grids of (AGN and star- 
formation) photoionization models now available should 
allow a new statistical approach to analyzing the contri- 
butions of shocks, star- formation and AGN photoioniza- 
tion to excitation of emission lines in galaxies. 

The complete set of electronic files that comprise the 
library are available on-line, along with tools to assist in 
the comparison of observations to the model predictions 
line ratio diagrams. Together this library provides one 
of the largest databases of radiative shock models and 
a unique tool in the interpretation and diagnosis of fast 
shocks. 
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TABLE 3 
Model Parameters 



Model name 


B (nG) 


B/VH (,uG 


cm 3 / 2 ) 




Solar abundance, n=1.0, 


v=100,125,.. 


.1000 


M_nl_bO 


1.0e-04 




1 On f\A 


M_nl_b0.5 


0.50 




n <^n 

U-OU 


M_nl_bl 


1.00 




i.UU 


M_nl.b2 


2.00 




o nn 


M_nl_be 


3.23 




o.zo 


M_nl_b4 


4.00 




4.00 


M_nl_b5 


5.00 




5.00 


M_nl_blO 


10.0 




10.0 


Solar x2 abundance, n=1.0, v=100,125 


,...1000 


R-nl.bO 


1.0e-04 




1.0e-04 


R.nl.b0.5 


0.50 




0.50 


R_nl_bl 


1.00 




1.00 


R_nl_b2 


2.00 




2.00 


R_nl_be 


3.23 




3.23 


R_nl_b4 


4.00 




4.00 


R.nl.b5 


5.00 




5.00 


R_nl_blO 


10.0 




10.0 


Dopita 2005 abundance, n= 


1.0, v=100,125,...1000 


J.nl.b0 


1.0e-04 




1.0c- 04 


J.nl_b0.5 


0.50 




0.50 


J_nl_bl 


1.00 




1.00 


J_nl_b2 


2.00 




2.00 


J_nl_be 


3.23 




3.23 


J_nl_b4 


4.00 




4.UU 


J_nl_b5 


5.00 




5.00 


J_nl_bl0 


10.0 
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Model name 


B (juG) 
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TABLE 5 

Integrated Column Densities for the precursor component of solar 
abundance models, n=1.0, b=3.23 



V=200 H He C N O Ne Mg Al Si S Ar Ca Fc 

I 8.836E+19 8.699E+18 4.030E+15 9.792E+15 7.582E+16 8.448E+15 1.350E+14 2.057E+11 1.271E+12 1.935E+11 1.637E+14 3.158E+12 1.500E+12 

II 1.198E+20 1.127E+19 4.675E+16 7.516E+15 7.372E+16 9.629E+15 4.418E+15 3.741E+14 5.617E+15 1.536E+15 2.199E+14 2.774E+13 2.343E+15 

III 0. 3.756E+17 2.471E+16 5.905E+15 2.685E+16 7.473E+15 3.305E+15 1.799E+14 1.247E+15 1.765E+15 3.684E+14 2.278E+14 1.594E+15 

IV 0. 0. 1.052E+14 1.471E+14 8.204E+14 6.430E+13 5.709E+13 5.960E+13 5.124E+14 7.515E+13 4.014E+12 1.520E+13 9.311E+14 

V 0. 0. 9.281E+11 4.191E+11 4.332E+12 2.582E+11 1.521E+11 6.680E+11 9.974E+12 4.006E+11 5.185E+09 5.757E+11 1.105E+13 

VI 0. 0. 1.429E+05 9.819E+08 3.015E+09 1.112E+08 1.138E+08 5.264E+08 1.494E+10 1.158E+10 1.928E+07 1.834E+09 6.478E+10 

VII 0. 0. 0. 0. 0. 0. 0. 2.240E+04 3.348E+05 8.522E+05 2.597E+04 6.880E+05 1.057E+08 



V=500 


H 


He 


c 


N 


o 


Ne 


Mg 


Al 


Si 


S 


Ar 


Ca 


Fe 




I 


3.923E+21 


3.676E+20 


9.843E+15 


4.300E+17 


3.322E+18 


2.010E+17 


3.353E+14 


3.900E+11 


2.609E+12 


4.559E+11 


6.127E+15 


8.244E+11 


1.177E+13 




II 


2.134E+21 


1.688E+20 


1.490E+18 


4.651E+16 


2.496E+17 


2.103E+17 


1.361E+17 


1.280E+16 


1.510E+17 


4.937E+16 


9.052E+15 


4.467E+13 


9.166E+16 




III 


0. 


5.563E+19 


3.977E+17 


8.208E+16 


1.088E+18 


2.747E+17 


5.874E+16 


1.473E+14 


1.658E+15 


2.348E+16 


2.249E+15 


6.091E+15 


8.385E+15 




IV 


0. 


0. 


1.280E+17 


9.142E+16 


3.732E+17 


4.007E+16 


1.562E+16 


1.884E+15 


1.431E+16 


1.683E+16 


3.971E+15 


1.158E+15 


3.125E+16 




V 


0. 


0. 


1.596E+17 


1.415E+16 


8.908E+16 


1.707E+16 


1.177E+16 


1.954E+15 


2.236E+16 


4.867E+15 


5.156E+14 


5.959E+14 


8.587E+15 




VI 


0. 


0. 


1.412E+16 


1.475E+16 


3.122E+16 


2.027E+15 


5.601E+15 


9.322E+14 


1.997E+16 


3.629E+15 


7.041E+13 


8.656E+13 


1.861E+15 




VII 


0. 


0. 


2.303E+14 


7.781E+14 


2.508E+15 


9.202E+12 


1.897E+15 


1.379E+14 


5.318E+15 


4.522E+13 


8.022E+12 


7.946E+12 


2.311E+14 




VIII 


0. 


0. 


0. 


3.292E+12 


2.088E+13 


2.861E+H 


2.077E+14 


2.284E+13 


3.325E+14 


1.316E+13 


9.636E+10 


5.233E+11 


1.628E+13 




IX 


0. 


0. 


0. 


0. 


1.810E+10 


6.851E+09 


4.352E+12 


8.903E+11 


1.887E+13 


8.379E+11 


7.297E+08 


1.572E+10 


6.101E+11 




X 


0. 


0. 


0. 


0. 


0. 


0. 


5.528E+10 


8.765E+09 


3.715E+11 


1.191E+10 


5.742E+07 


8.657E+07 


2.423E+10 




XI 


0. 


0. 


0. 


0. 


0. 


0. 


2.444E+08 


4.969E+07 


1.768E+09 


2.439E+08 


1.302E+06 


1.864E+05 


9.224E+07 




XII 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


1.431E+05 


5.939E+06 


2.380E+06 


8.874E+03 


4.899E+03 


4.744E+04 




XIII 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


1.078E+04 


5.912E+03 


69.6 


23.8 


78.0 


> 


XIV 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


11.5 


0.132 


2.108E-02 


9.953E-02 


XV 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


2.409E-05 


0. 


ST 


V=1000 


H 


He 


c 


N 


o 


Ne 


Mg 


Al 


Si 


S 


Ar 


Ca 


Fc 




I 


3.711E+22 


3.599E+21 


5.026E+15 


4.214E+18 


3.187E+19 


2.568E+18 


1.684E+14 


2.034E+11 


1.355E+12 


1.640E+11 


7.177E+16 


5.669E+09 


6.059E+12 


CD 
r+ 


II 


1.165E+22 


6.939E+20 


1.415E+19 


2.041E+17 


1.591E+18 


1.583E+18 


1.272E+18 


1.167E+17 


1.394E+18 


5.163E+17 


7.005E+16 


9.875E+12 


8.638E+17 




III 


0. 


4.722E+20 


6.107E+17 


8.576E+16 


2.672E+18 


1.105E+18 


2.865E+17 


1.564E+14 


2.784E+15 


1.231E+17 


3.955E+15 


5.201E+16 


5.999E+16 




IV 


0. 


0. 


3.272E+17 


2.090E+17 


1.041E+18 


1.308E+17 


4.236E+16 


3.643E+15 


1.659E+16 


4.578E+16 


1.681E+16 


4.763E+15 


1.064E+17 




V 


0. 


0. 


1.055E+18 


8.139E+16 


1.134E+18 


2.714E+17 


3.561E+16 


5.219E+15 


5.125E+16 


1.984E+16 


2.512E+15 


1.945E+15 


1.720E+16 




VI 


0. 


0. 


1.284E+18 


3.894E+17 


1.525E+18 


3.198E+17 


4.021E+16 


4.439E+15 


7.136E+16 


2.582E+16 


3.146E+15 


1.468E+15 


1.826E+16 




VII 


0. 


0. 


2.683E+17 


2.619E+17 


1.260E+18 


1.640E+16 


7.340E+16 


3.461E+15 


7.833E+16 


4.853E+15 


7.444E+15 


1.755E+15 


3.025E+16 




VIII 


0. 


0. 


0. 


2.588E+16 


3.918E+17 


4.008E+15 


7.667E+16 


5.970E+15 


5.138E+16 


2.269E+16 


9.964E+14 


1.680E+15 


2.729E+16 




IX 


0. 


0. 


0. 


0. 


1.916E+16 


1.325E+15 


2.119E+16 


3.668E+15 


4.620E+16 


2.376E+16 


1.095E+14 


6.088E+14 


1.274E+16 




X 


0. 


0. 


0. 


0. 


0. 


3.607E+13 


5.216E+15 


6.005E+14 


1.639E+16 


6.145E+15 


1.663E+14 


3.904E+13 


6.620E+15 




XI 


0. 


0. 


0. 


0. 


0. 


1.792E+11 


4.759E+14 


7.096E+13 


1.380E+15 


2.123E+15 


7.028E+13 


1.185E+12 


4.243E+14 




XII 


0. 


0. 


0. 


0. 


0. 


0. 


2.139E+12 


3.989E+12 


9.250E+13 


3.101E+14 


8.348E+12 


8.244E+11 


4.262E+12 




XIII 


0. 


0. 


0. 


0. 


0. 


0. 


1.609E+09 


4.506E+09 


3.189E+12 


1.289E+13 


1.452E+12 


1.142E+11 


1.425E+11 




XIV 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


1.494E+09 


6.780E+11 


8.089E+10 


4.113E+09 


4.040E+09 




XV 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


1.078E+10 


1.054E+09 


1.893E+08 


6.014E+07 




XVI 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


1.465E+07 


3.442E+06 


1.810E+04 




XVII 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


7.387E+04 


1.705E+04 


28.9 




XVIII 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


58.0 


0.290 




XIX 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


5.533E-02 


7.180E-04 




XX 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


5.351E-07 




XXI 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


5.694E-10 





The MAPPINGS III Shock Model Library 



37 



TABLE 6 

Line Ratios for shock components of solar abundance models, n=1.0, 

B=3.23 
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2.058 




O n] AA3727, 3729 


10.652 


5.964 


11.045 


14.156 


14.601 


14.141 


14.230 


14.770 


15.428 


15.446 




S n] AA4067, 4076 


0.066 


0.057 


0.109 


0.138 


0.133 


0.144 


0.160 


0.177 


0.189 


0.200 




O ml A4363 


0.290 


0.115 


0.073 


0.061 


0.055 


0.050 


0.045 


0.044 


0.045 


0.045 


He ii A4686 


0.039 


0.060 


0.154 


0.452 


0.337 


0.269 


0.249 


0.253 


0.273 


0.274 


H 8 A4861 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 




O m] A5007 


3.190 


1.188 


1.253 


2.688 


3.045 


3.060 


3.039 


3.167 


3.490 


3.714 




N i] A5200 


0.109 


0.080 


0.242 


0.292 


0.462 


0.968 


1.428 


1.793 


2.005 


2.206 




Fc vn] A6085 


0.000 


0.013 


0.015 


0.010 


0.008 


0.006 


0.005 


0.004 


0.004 


0.003 




O L] A6300 


0.145 


0.242 


0.896 


1.225 


1.531 


2.467 


3.365 


4.112 


4.599 


4.998 


H a A6563 


3.214 


3.010 


2.938 


2.936 


2.947 


2.968 


2.999 


3.029 


3.060 


3.082 




N n] A6583 


1.641 


2.024 


4.155 


5.103 


5.056 


4.871 


4.870 


5.040 


5.265 


5.269 




S n] AA6716, 6731 


1.445 


1.975 


2.919 


3.200 


3.050 


3.504 


4.066 


4.543 


4.801 


5.074 




O n] AA7318, 7324 


0.449 


0.170 


0.195 


0.247 


0.264 


0.261 


0.272 


0.292 


0.318 


0.326 




S ml AA9069, 9532 


0.751 


0.840 


2.421 


3.227 


3.896 


4.871 


5.672 


6.311 


6.763 


7.046 


Br a A4.051^m 


0.112 


0.096 


0.087 


0.084 


0.085 


0.086 


0.087 


0.087 


0.086 


0.086 




Ar vi] A4.530^m 


0.000 


0.000 


0.000 


0.000 


9.454e-5 


7.785e-5 


6.537e-5 


5.679c-5 


A n A o „ er 

4.948C-5 


4.345c-5 




Ar ii] A6.983^m 


0.089 


0.092 


0.080 


0.080 


0.107 


0.171 


0.222 


0.252 


0.267 


0.280 




Ar in] A6.983^tm 


0.038 


0.114 


0.300 


0.339 


0.334 


0.302 


0.273 


0.256 


0.248 


0.239 




Nc vi] A7.652^m 


0.000 


0.014 


0.011 


0.007 


0.005 


0.004 


0.004 


0.003 


0.003 


0.002 




S iv] A10.44/un 


0.022 


0.011 


0.016 


0.041 


0.053 


0.064 


0.077 


0.093 


0.112 


0.129 




Ne ii] A12.8/im 


0.628 


0.840 


0.741 


0.569 


0.907 


1.543 


2.028 


2.224 


2.293 


2.390 




Ne v] A14.5^tm 


0.000 


0.013 


0.007 


0.004 


0.004 


0.003 


0.003 


0.002 


0.002 


0.002 




Ne III] A15.5/mi 


0.310 


0.446 


2.230 


3.229 


3.435 


3.936 


4.427 


4.799 


5.088 


5.291 




S m] A18.7Mm 


0.269 


0.440 


1.052 


1.298 


1.647 


2.094 


2.402 


2.594 


2.712 


2.783 




O iv] A25.9/jm 


0.083 


0.134 


0.097 


0.347 


0.344 


0.348 


0.370 


0.433 


0.556 


0.655 


logio H/3 (erg cm * 


s- 1 ) -5.317 


-4.444 


-3.996 


-3.675 


-3.469 


-3.303 


-3.155 


-3.032 


-2.918 


-2.814 
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TABLE 7 

Line Ratios for precursor components of solar abundance models, 

n=1.0, B=3.23 



Line Shock Velocity (km s 1 ) 







200 


300 


400 


500 


600 


700 


800 


900 


1000 


C in A977 


0.000 


7.075c-5 


0.001 


0.009 


0.075 


0.200 


0.275 


0.288 


0.285 


N in A991 


0.000 


1.013e-5 


0.000 


0.002 


0.015 


0.038 


0.053 


0.050 


0.050 


O vi AA1032, 1037 


0.000 


0.000 


0.000 


0.011 


0.653 


4.650 


12.084 


22.156 


30.910 


Ly a A1215 


19.540 


19.610 


20.190 


22.220 


26.280 


30.870 


34.090 


37.440 


41.390 


N v A1240 


0.000 


0.000 


6.491e-5 


0.014 


0.217 


0.673 


1.055 


1.333 


1.468 


Si iv + O iv A1400 


1.099c-5 


0.002 


0.012 


0.117 


0.689 


1.396 


1.831 


2.089 


2.464 


N iv A1486 


0.000 


0.000 


0.002 


0.047 


0.241 


0.431 


0.515 


0.547 


0.570 


C iv A1550 


0.000 


0.003 


0.035 


1.206 


5.764 


10.272 


12.136 


12.485 


12.526 


He II A1640 


0.263 


1.425 


3.617 


2.246 


2.091 


2.358 


2.625 


3.033 


3.162 


O in] A1664 


9.984e-5 


0.005 


0.030 


0.150 


0.558 


1.040 


1.430 


1.638 


1.892 


N in A1750 


5.054e-5 


0.001 


0.006 


0.021 


0.054 


0.076 


0.080 


0.065 


0.058 


C in] A1909 


0.004 


0.085 


0.298 


1.022 


2.113 


2.806 


2.860 


2.630 


2.408 


C n] A2326 


0.012 


0.036 


0.063 


0.137 


0.210 


0.266 


0.293 


0.324 


0.381 


Mg II A2800 


0.107 


0.290 


0.410 


0.716 


0.923 


1.020 


1.032 


1.047 


1.111 




Nc v] A3426 


0.000 


0.001 


0.004 


0.147 


0.686 


1.484 


2.089 


2.631 


2.825 




Ne m] A3869 


0.052 


0.310 


0.569 


1.310 


2.070 


2.471 


2.636 


2.584 


2.569 




O n] AA3727, 3729 


1.119 


1.241 


1.065 


1.160 


1.316 


1.416 


1.415 


1.416 


1.526 




S o] AA4067, 4076 


0.005 


0.006 


0.007 


0.014 


0.018 


0.022 


0.026 


0.030 


0.036 




ml A4363 


0.000 


0.009 


0.034 


0.120 


0.298 


0.448 


0.549 


0.586 


0.638 


He II A4686 


0.045 


0.234 


0.603 


0.340 


0.300 


0.327 


0.357 


0.407 


0.420 


H /3 A4861 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


[0 in] A5007 


0.661 


4.594 


7.639 


16.270 


23.960 


26.920 


27.860 


26.760 


26.760 


[N i] A5200 


0.004 


0.021 


0.057 


0.142 


0.203 


0.255 


0.294 


0.343 


0.420 


[Fe vn] A6085 


0.000 


0.000 


7.858e-5 


0.003 


0.021 


0.078 


0.142 


0.199 


0.205 


[O i] A6300 


0.013 


0.062 


0.131 


0.316 


0.463 


0.571 


0.633 


0.708 


0.829 


H a A6563 


3.011 


2.961 


2.943 


2.873 


2.876 


2.899 


2.922 


2.947 


2.976 


[N n] A6583 


0.650 


0.465 


0.325 


0.375 


0.442 


0.475 


0.460 


0.443 


0.449 


[S n] AA6716, 6731 


0.244 


0.165 


0.180 


0.334 


0.426 


0.526 


0.613 


0.715 


0.855 


[0 n] AA7318, 7324 


0.006 


0.014 


0.016 


0.019 


0.024 


0.027 


0.027 


0.027 


0.029 


[S in] AA9069, 9532 




1 9QQ 


u. yuo 


1 014 




l.UUO 


u.you 




922 


Br a A4.051um 


0.102 


0.094 


0.092 


0.081 


0.076 


0.073 


0.072 


0.070 


0.069 




Ar vi] A4.530/im 


0.000 


3.302e-5 


0.001 


0.003 


0.008 


0.019 


0.028 


0.033 


0.034 




Ar n] A6.983/im 


0.013 


0.006 


0.008 


0.015 


0.018 


0.022 


0.024 


0.027 


0.031 




Ar in] A6.983/im 


0.264 


0.282 


0.159 


0.130 


0.109 


0.090 


0.075 


0.060 


0.051 




Ne vi] A7.652^m 


0.000 


0.000 


0.004 


0.031 


0.138 


0.373 


0.633 


0.929 


1.153 




S iv] A10.44/mi 


0.144 


1.993 


3.005 


3.303 


3.149 


2.851 


2.588 


2.298 


2.167 




Ne n] A12.8^m 


0.357 


0.091 


0.066 


0.107 


0.116 


0.133 


0.155 


0.181 


0.216 




Nc v] A14.5/im 


1.702e-5 


0.007 


0.036 


0.098 


0.196 


0.302 


0.363 


0.417 


0.419 




Ne III] A15.5^tm 


1.147 


2.148 


2.006 


2.321 


2.390 


2.281 


2.143 


1.947 


1.849 




S in] A18.7/mi 


0.848 


0.667 


0.367 


0.334 


0.315 


0.300 


0.284 


0.273 


0.283 




O iv] A25.9Atm 


0.342 


5.524 


15.890 


12.550 


12.190 


11.330 


10.150 


9.131 


8.604 


logio H/3 (erg cm ' 2 s x ) 


-4.622 


-4.130 


-3.769 


-3.599 


-3.464 


-3.359 


-3.261 


-3.164 


-3.082 
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TABLE 8 

Line Ratios for shock + precursor components of solar abundance 

MODELS, N=1.0, B = 3.23 





Line 








Shock Velocity (k 


n s j 












zuu 


ouu 


^tuu 


500 


600 


1 uu 


ouu 


yuu 


1UUU 


C in A977 


0.593 


0.267 


0.174 


0.144 


0.146 


0.176 


0.189 


0.182 


0.169 


N in A991 


0.373 


0.201 


0.125 


0.101 


0.089 


0.086 


0.082 


0.072 


0.066 


O VI AA1032, 1037 


30.700 


16.634 


12.392 


10.507 


8.917 


9.088 


10.818 


13.538 


15.709 


Ly a A1215 


25.055 


24.923 


27.313 


29.194 


32.916 


38.213 


43.695 


49.091 


53.477 


N v A1240 


0.810 


0.448 


0.289 


0.225 


0.265 


0.108 


0.521 


0.596 


0.615 


Si iv + O iv A1400 


1.448 


0.750 


0.517 


0.454 


0.621 


0.830 


0.941 


0.993 


1.079 


N iv A 1486 


0.073 


0.037 


0.028 


0.040 


0.115 


0.180 


0.203 


0.209 


0.210 


C iv A1550 


2.759 


1.306 


0.791 


1.113 


2.843 


4.366 


4.866 


4.841 


4.674 


He ii A1640 


0.393 


1.186 


3.231 


2.222 


1.904 


1.920 


2.029 


2.254 


2.292 


O in] A1664 


0.427 


0.247 


0.182 


0.208 


0.354 


0.511 


0.632 


0.688 


0.752 


N in A1750 


0.051 


0.030 


0.024 


0.027 


0.038 


0.044 


0.044 


0.037 


0.034 


C in] A1909 


0.436 


0.367 


0.586 


0.925 


1.345 


1.583 


1.612 


1.581 


1.511 


C n] A2326 


0.268 


0.338 


0.414 


0.492 


0.650 


0.845 


1.029 


1.176 


1.301 


Mg ii A2800 


0.381 


0.898 


1.149 


1.426 


1.851 


2.193 


2.561 


2.832 


3.199 




Nc v] A3426 


0.173 


0.091 


0.056 


0.107 


0.318 


0.604 


0.805 


0.981 


1.015 




Ne in] A3869 


0.207 


0.532 


0.948 


1.311 


1.682 


1.921 


2.089 


2.184 


2.237 




O n] AA3727, 3729 


4.030 


6.897 


8.321 


8.882 


8.902 


9.304 


9.815 


10.352 


10.567 




S n] AA4067, 4076 


0.036 


0.066 


0.080 


0.082 


0.092 


0.107 


0.121 


0.132 


0.142 




O ml A4363 


0.070 


0.046 


0.049 


0.083 


0.151 


0.200 


0.231 


0.241 


0.253 


He ii A4686 


0.054 


0.188 


0.519 


0.339 


0.282 


0.279 


0.292 


0.322 


0.325 


H /3 A4861 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


[0 in] A5007 


0.978 


2.667 


4.895 


8.672 


11.598 


12.220 


12.328 


11.920 


11.792 


[N i] A5200 


0.050 


0.148 


0.187 


0.326 


0.655 


0.977 


1.237 


1.403 


1.580 


[Fe vn] A6085 


0.008 


0.009 


0.005 


0.005 


0.012 


0.033 


0.056 


0.075 


0.074 


[O i] A6300 


0.151 


0.543 


0.737 


1.014 


1.648 


2.291 


2.821 


3.189 


3.537 


H a A6563 


3.010 


2.948 


2.939 


2.916 


2.930 


2.961 


2.989 


3.019 


3.045 


[N n] A6583 


1.476 


2.593 


2.973 


3.064 


3.062 


3.180 


3.341 


3.518 


3.580 


[S n] AA6716, 6731 


1.284 


1.754 


1.854 


1.895 


2.246 


2.705 


3.085 


3.321 


3.595 


[O n] AA7318, 7324 


0.104 


0.118 


0.144 


0.160 


0.164 


0.177 


0.194 


0.213 


0.222 


[S m] AA9069, 9532 


0.926 


1.946 


2.219 


2.670 


3.300 


3.877 


4.322 


4.640 


4.899 


Br a A4.051/im 


0.099 


0.090 


0.088 


0.083 


0.082 


0.082 


0.081 


0.080 


0.080 




Ar vi] A4.530Mm 


0.000 


0.000 


0.001 


0.001 


0.003 


0.008 


0.011 


0.012 


0.012 




Ar ii] A6.983^tm 


0.060 


0.049 


0.048 


0.068 


0.108 


0.145 


0.167 


0.180 


0.192 




Ar in] A6.983Mm 


0.174 


0.292 


0.259 


0.247 


0.223 


0.203 


0.189 


0.179 


0.173 




Ne vi] A7.652^tm 


0.008 


0.007 


0.005 


0.016 


0.059 


0.146 


0.237 


0.338 


0.406 




S iv] A10.44^m 


0.064 


0.852 


1.362 


1.436 


1.324 


1.143 


1.018 


0.904 


0.844 




Nc ii] A12.8^m 


0.647 


0.466 


0.345 


0.567 


0.960 


1.299 


1.457 


1.528 


1.628 




Ne v] A14.5/im 


0.008 


0.007 


0.018 


0.044 


0.082 


0.117 


0.136 


0.152 


0.148 




Ne in] A15.5/im 


0.726 


2.195 


2.684 


2.961 


3.304 


3.602 


3.814 


3.950 


4.085 




S m] A18.7^m 


0.603 


0.889 


0.883 


1.088 


1.367 


1.594 


1.737 


1.828 


1.907 




O iv] A25.9^m 


0.217 


2.393 


7.275 


5.538 


5.186 


4.584 


4.038 


3.662 


3.441 



logi H/3 (erg cm" 2 s" 1 ) -4.223 -3.757 -3.418 -3.228 -3.075 -2.944 -2.831 -2.723 -2.627 



